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|SS  ABSTRACT  I 

A  uniformly  heated,  pressure-driven  shock  tube  was  used  to  produce  4200  K  to  12,000  K 
mercury  plasmas  in  which  the  relaxation  process  was  studied.  The  total  plasma 
luminosity  was  observed  using,  a  photomultiplier  tube,  with  the  time  of  shock  passage 
and  shock  velocity  establishes  using  a  schlieren  optical  system.  It  is  shown  that  the 
ionizational  relaxation  time  depends  inversely  on  the  initial  vapor  pressure  or  atom 
number  density,  and  directly  on  the  final  electron  number  density.  The  final  plasma 
temperature  is  shown  to  be  more  relevant  for  strongly  ionized  shocks  than  the  commonly 
used  ideal  gas  temperature.  Ionizations  caused  by  electron-atom  collisions  appear  to 
occur  as  a  single  step  rather  than  by  cumulative  steps  as  atom-atom  ionizations.  A 
previously  unobserved  feature  of  ionizational  relaxation  has  been  detected.  The 
luminosity  profile  shows  two  distinct  points  at  which  the  slope  of  the  curve  changes 
abruptly.  No  explanation  of  this  feature  could  be  found.  Radiation  overshoot  is 
shown  to  be  an  effect  which  occurs  at  low  impurity  levels.  The  true  relaxation  time 
corresponds  to  the  end  of  the  overshoot  phenomenon.  Impurities  are  seen  to  affect 
the  shape  of  the  relaxation  profile  greatly,  and  manometer  and  vacuum  pump  oils  are 
shown  to  be  possible  sources  of  impurity  in  some  previous  studies. 
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ABSTRACT 


This  paper  is  concerned  with  aspects  of  the  ionizational 
relaxation  in  a  shock-heated  mercury  plasma.  A  uniformly 
heated,  pressure-driven  shock  tube  was  used  to  produce 
4200  °K  to  12,000  °K  mercury  plasmas  in  which  the  relaxation 
process  was  studied.  The  total  plasma  luminosity  was  ob¬ 
served  using  a  photomultiplier  tube,  with  the  time  of  shock 
passage  and  the  shock  velocity  established  using  a  schlieren 
optical  system. 

It  is  shown  that  the  ionizational  relaxation  time 
depends  inversely  on  the  initial  vapor  pressure  or  atom 
number  density,  and  directly  on  the  final  electron  number 
density.  The  f^nal  plasma  temperature  is  shown  to  be  more 
relevant  for  strongly  ionized  shocks  than  the  commonly  used 
ideal  gas  temperature. 

Ionizations  caused  by  electron-atom  collisions  appear 
to  occur  as  a  single  step  rather  than  by  cumulative  steps 
as  atom-atom  collision  ionizations. 

A  previously  unobserved  feature  of  ionizational  relax¬ 
ation  has  been  detected.  The  luminosity  profile  shows  two 
distinct  points  at  which  the  slope  of  the  curve  changes 
abruptly.  No  explanation  of  this  feature  could  be  found. 

Radiation  overshoot  is  shown  to  be  an  effect  which 
occurs  at  low  impurity  levels.  The  true  relaxation  time 
corresponds  to  the  end  of  the  overshoot  phenomenon.  Im¬ 
purities  are  seen  to  affect  the  shape  of  the  relaxation  pro¬ 
file  greatly,  and  manometer  and  vacuum  pump  oils  are  shown 
to  be  possible  sources  of  impurities  in  some  previous  studies. 
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CHAPTER  I 


INTRODUCTION 

1.  Elementary  Shock  Theory 

A  simple  pressure-driven  shock  tube  consists  of  a  long 
tube  of  uniform  cross-section,  divided  into  two  chambers 
by  a  breakable  diaphragm.  The  expansion  chamber  contains  a 
test  gas  at  low  pressure  (1  to  100  Torr  usually).  The 
compression  chamber  is  slowly  filled  with  a  "driver"  gas 
until  the  diaphragm  bursts  (100  to  2000  psi)  and  the  driver 
Is  released  into  the  expansion  chamber-.  Diaphragm  thickness, 
material,  and  the  depth  of  grooves  scored  on  the  diaphragm 
face  control  the  bursting  pressure. 

If  the  expansion  of  the  driver  into  the  low  pressure 
chamber  is  supersonic  relative  to  the  i-est  gas,  then  the  gas 
far  ahead  of  the  driver  remains  at  rest  in  its  initial  state 
(P0,  T0,yQj)i  until  the  instant  it  is  "shocked".  Diffusion 
between  the  two  gases  is  so  extremely  slow  compared  to  the 
shock  processes  that  all  mixing  ii  *y  normally  be  ignored. 
Therefore,  the  interface  between  the  gases  acts  similar  to 
a  solid  piston,  moving  at  uniform  speed  and  sweeping  up  the 
test  gas  ahead  of  it. 

The  test  gas  is  compressed  and  accelerated  until  it  has 
the  same  pressure  and  speed  as  the  interface,  but  the  actual 
pressure  change  in  the  test  gas  moves  forward  even  faster 
than  the  interface.  This  sudden  pressure  step  in  the  test 
gas  is  the  shock  front. 
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When  the  diaphragm  hursts,  the  desired  single  sharp 
step  in  pressure,  extending  uniformly  across  the  entire  cross- 
section,  will  not  normally  occur.  A  single  plane  pressure 
discontinuity  will  result,  however,  after  the  pressure 
change  has  moved  forward  20  to  pO  tube  diameters.  Pressure 
variations  perpendicular  to  the  flow  will  obviouslv  come 
to  an  equilibrium  very  quickly  as  the  interface  and  shock 
front  move  forward.  Pressure  variations  along  the  tube  axis 
will  act  as  a  series  of  pressure  pulses  following  one 
another.  The  gas  will  be  heated  by  the  rapid  compression 
due  to  the  passage  of  the  first  pressure  pulse.  The  vel¬ 
ocity  of  sound  (or  pressure  change)  increases  as  the  square 
root  of  the  absolute  temperature,  so  that,  a  following 
pulse  will  travel  faster  than  the  first.  Therefore,  a 
second  pulse  will  overtake  the  first  and  combine  with  it 
to  form  a  single  large  pulse.  Similarly,  all  gradual 
variations  in  pressure  will  steepen  as  the  tail  of  the 
variation  overtakes  the  head,  until  a  pressure  discontin¬ 
uity  results.  A  more  detailed  description  of  how  irregular¬ 
ities  in  pressure  converge  to  form  the  shock  front  has  ceen 

i 

given  by  Laporte  . 

In  the  compression  section  s  rarefaction  wave  travels 
back  into  the  driver  gas,  but  because  conditions  here  are 
the  reverse  of  those  occurring  in  the  test  section,  this 
wave  gradually  disperses  along  the  tube  as  it  moves. 

The  "shocked"  gas,  which  has  all  undergone  exactly  the 
same  compression  process  in  reaching  its  final  state 
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(P,  ,  T,  ,p(  ),  forms  a  slug  of  gas  at  uniform  high  temper¬ 
ature  and  pressure.  This  slug  of  plasma,  between  the  ■‘.iter- 
face  and  the  shock  front,  moves  at  the  speed  of  the  inter¬ 
face  and  lengthens  as  more  gas  collects  at  the  shock  front. 
Because  there  Is  a  Boltzmann  distribution  over  the  excited 
atomic  and  molecular  states,  and  an  equilibrium  of  energy 
among  the  degrees  of  freedom  of  the  system,  this  plasma  is 
usually  referred  to  as  being  in  "local  thermodynamic 
equilibrium"  (L.T.E. ). 

The  uniformity  of  this  plasma  and  the  existence  of  L.T.E. 
make  it  excellent  as  a  spectroscopic  source,  or  for  studies  of 

O 

supersonic  flows  .  Measurements  of  transition  probabilities 

(gf  values)  for  optical  frequencies^,  and  the  broadening 
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of  spectral  lines  *  ,  have  been  particularly  successful 
using  shock  tubes. 

To  the  particles  in  the  plasma,  the  tube  end-wall  appears 
as  a  piston  rushing  toward  them  at  supersonic  speed.  There¬ 
fore,  when  this  pl«uma  slug  meets  the  end-wall  a  new  shock 
front  will  form  and  travel  back  thr  ugh  the  slug,  compress¬ 
ing  and  heating  the  gas  still  further.  Similarly,  this 
shock  wave  w^ll  be  reflected  when  it  strikes  the  interface. 
This  second  reflection  is  not  generally  as  good  as  the  first, 
however,  because  the  interface  is  slightly  diffuse  due  to 
some  mixing  of  the  driver  and  test  gases. 

The  gas  plasma  fin  L.T.E.)  behind  each  shock  may  be  used 
for  study  purposes,  but  the  time  for  observation  gets  pro¬ 
gressively  shorter  behind  each  reflection.  The  observation 


*> 

time  available  is  from  the  passage  of  the  shock  front  until 
the  arrival  of  either  the  interface  or  the  next  reflected 
shock.  For  this  reason,  studies  are  commonly  limited  to 
the  primary  or  first  reflected  shock  regions. 

Figure  1  illustrates  all  of  these  actions  graphically  . 

2.  Boundary  Layers 

For  the  hot  plasma  moving  at  supersonic  speed,  energy 
loss  by  the  gas  molecules  and  viscous  drag  occur  at  the 
cold  tube  wall.  The  combination  of  cooling  and  drag  results 
in  the  growth  of  a  cooler,  slower  movirg  "boundary  layer" 
along  each  wall,  and  shock  front  deceleration  is  possible. 

The  boundary  flow  is  usually  very  thin  and  laminar 
near  the  shock  front,  and  becomes  thicker  and  turbulent 
nearer  to  the  interface.  As  shock  strength  increases  the 
layer  will  often  become  thicker,  and  turbulence  will  begin 
earlier.  Boundary  layers  in  mercury  vapor  have  been  ob¬ 
served  to  grow  to  such  an  extent  that  the  tube  is  complete¬ 
ly  choked  off  and  only  a  supersonic  turbulent  flow  results. 

The  available  observation  time  will  be  shortened  because 
material  is  removed  from  the  plasma  and  flows  back  through 
the  boundary  layer  to  mix  with  the  driver  gas.  When  ob¬ 
servations  are  analyzed,  the  presence  of  a  boundary  layer 
must  not  be  forgotten,  because  observed  phenomena  may 
actually  be  in  this  layer,  not  in  the  uniform  plasma  being 


studied. 
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3.  Radiation 

For  extremely  luminous  shocks  the  energy  lost  by  radi¬ 
ation  can  become  very  significant  and  true  thermodynamic 
equilibrium  may  not  actually  occur.  Since  the  plasma  is 
optically  quite  thin,  radiation  easily  escapes  from  the 
system  with  a  resulting  rapid  energy  loss.  The  temperature 
and  luminosity  decrease  toward  the  tail  of  the  slug,  accom¬ 
panied  by  adjustments  in  the  other  gas  variables  such  as 
density  and  number  of  free  electrons,  etc.  The  decrease  in 
the  total  energy  of  the  plasma  gas  with  time,  caused  by 
radiation  losses  is  known  as  radiation  decay. 

Radiation  decay  has  been  observed  for  many  gases,  either 

for  very  strong  primary  shocks ,  or  for  the  reflected  shocks. 

7 

Tumakaev  and  Lazovskaya  have  studied  the  very  pronounced 
decay  which  occurs  fo"  only  moderately  strong  shocks  in 
mercury  vapor.  Their  observations  will  be  discussed  further 
in  connection  with  the  present  study  of  mercury  vapor. 

4.  Relaxations 

Ideally  the  passage  of  the  shock  front  Instantaneously 
adds  energy  and  compresses  the  stationary  test  gas,  causing 

the  gas  to  change  from  its  initial  condition  ( Pe  ,  p0  ,  T0  ) 
to  some  new  equilibrium  condition  However, 

the  only  efficient  mechanism  available  to  distribute  the  add¬ 
ed  energy  among  the  available  energy  states  and  bring  about 
final  equilibrium,  is  collisions  between  gas  molecules. 

The  time  between  collisions  is  not  infinitesimal,  and 
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many  collisions  will  be  necessary  to  bring  about  equilibrium, 
therefore,  the  change  ir.  the  state  of  the  gas  cannot  occur 
instantly.  This  delay  was  first  observed  at  Michigan  by 
Turner  in  1956^ ,  but  the  relaxation  processes  are  still  not 
fully  understood. 

Different  types  of  energy  states  (rotation,  vibration, 

O 

etc.)  actually  come  to  "equilibrium"  at  different  times  . 

The  "relaxation  time"  (T  )  is  defined  as  the  time  from  the 
passage  of  the  shock  front  until  equilibrium  for  that  type 
of  energy  state  occurs.  After  approximately  twice  the  mean 
time  between  collisions  a  Maxwellian  distribution  of  molecu¬ 
lar  velocities  will  exist.  This  is  the  relaxation  time  for 
■equilibrium  among  the  translational  energy  states. 

The  rotational  energy  states,  if  any,  are  found  to 
come  into  equilibrium  with  the  translational  states  after 
10  to  300  collisions  per  molecule,  depending  upon  the  gas. 
Next  the  vibrational  states  come  to  equilibrium  with  the 
rotational  and  translational  states.  Molecular  dissociation 
is  the  limit  cf  vibration  and  will  reach  equilibrium  at 
approximately  the  same  time.  After  a  few  thousand  collisions 
the  electronic  excitation  states  come  to  equilibrium  with 
the  others.  Ionization  represents  the  limiting  case  of 
electronic  excitation  and  occurs  almost  simultaneously. 

Ionization  relaxation  is  the  longest  and  therefore  the 
most  easily  studied  for  any  gas.  To  study  any  relaxation 
time,  the  gas  and  the  conditions  can  often  be  chosen  to 
suppress  or  lengthen  the  other  relaxations.  For  example, 
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rotational  relaxation  studies  have  been  done  using  diatomic 

9 

molecules  in  a  cooled  tube  with  weak  shocks  so  that  essen¬ 
tially  only  two  rotation  states  are  excited. 

During  the  relaxation  processes  the  pressure  will  remain 
approximately  constant  but  the  other  gas  variables  will  ad¬ 
just  toward  their  final  values.  The  particle  number  density 
or  mass  density  will  always  increase,  while  the  translational 
temperature  decreases.  Note  that  during  the  relaxation, 
"temperature"  has  no  clear  meaning  since  each  energy  type 
will  show  a  different  apparent  temperature  (or  temperature 
cannot  be  defined). 

5.  The  Problem 

In  1957  Petschek  and  Byron-10  published  the  results  of  a 
study  of  ionization  relaxation  times  in  argon  using  an 
electron  probe  to  measure  the  Instantaneous  density  of  elect¬ 
rons  in  the  plasma.  Their  conclusions  have  been  generally 
accepted  since  that  time,  but  several  questions  concerning 
their  treatment  of  results  will  be  discussed  in  Chapter  V. 

To  find  the  activation  energy  for  the  relaxation  process 
the  logarithm  of  the  relaxation  time  (log  7^  )  was  plotted 
versus  the  reciprocal  of  the  absolute  plasma  temperature 
(l/T).  A  parallel  straight  line  was  obtained  for  each 
initial  pressure,  P„  ,  showing  T  is  proportional  to 
exp  (E*/kT)  as  expected.  All  data  reduced  reasonably  well 
to  one  line  if  log  (P ,T  )  vs  (l/T)  was  plotted  instead. 

This  seemed  reasonable  since  collision  frequency  was  ex- 
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pected  to  be  proportional  to  the  pressure  and  'jf  to  be  in¬ 
versely  proportional  to  the  collision  frequency.  This 
experiment  is  accepted  as  evidence  that  the  ionization 
relaxation  time  is  inversely  proportional  to  the  pressure, 
or  number  density. 

In  1968,  as  part  of  a  study  of  mercury  vapor  in  a 

heated  shock  tube  Y.M.  Kim  included  a  brief  study  of  ion- 

6 

ization  relaxation0.  He  found  that  log  )  vs  (1/T)  was 

not  nearly  as  good  a  line  as  log  vs  (1/T).  It  therefore 
appeared  that  no  pressure  dependence  of  exists  for  mer¬ 
cury  vapor--in  direct  contradiction  to  the  results  for  argon. 
Recently  Dunaev  et  al"1-1  have  tried  to  account  for  this 
difference  by  postulating  the  formation  of  an  excited  di¬ 
atomic  mercury  molecule  which  acts  as  the  rate-controlling 
species  during  relaxation. 

Kim's  result  might  be  questioned  because  it  was  based 
on  only  13  points,  P#  was  varied  less  than  a  factor  of  2, 
and  the  shocks  at  the  small  T  end  of  the  graph  were  subject 
to  relatively  large  error.  However,  the  difference  in  pres¬ 
sure  dependence  may  be  due  to  other  inherent  differences  in 
the  gases  or  techniques.  Both  gases  are  monatomic  (above 
room  temperature)  but  the  low  ionization  potential  of  mer¬ 
cury  results  in  a  very  much  higher  final  degree  of  ion¬ 
ization.  Also,  the  massive  atoms  cause  shocks  in  mercury 
to  be  much  stronger  than  "similar  shocks"  in  argon  because 
shock  strength  depends  on  the  Mach  number  not  on  the  absolute 
velocity  (see  Chapter  II). 
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For  these  reasons.  It  was  felt  that  Kim's  conclusion 
should  be  checked  and  that  it  was  especially  important  to 
observe  much  weaker  shocks  in  mercury  as  a  better  comparison 
with  argon.  That  study  is  the  main  purpose  of  the  present 
work.  Also,  more  data  might  help  in  explaining  this  unex¬ 
pected  result  for  mercury  and  lead  to  a  better  understanding 
of  the  ionization  relaxation  process  as  well. 

It  has  also  been  apparent  for  some  time  that  minute 
amounts  of  impurities  have  a  very  major  effect  on  the  re¬ 
laxation  times.  This  was  noted  first  by  Petschek  and 
10 

Byron  ,  but  the  role  which  impurities  play  is  still  not 
understood.  One  great  advantage  of  the  heated  shock  tube 
is  the  ability  to  bake  out  the  tube  and  obtain  unusually 
high  test  gas  purity.  Therefore,  a  brief  study  of  impurity 
effects  is  included  in  this  present  work. 


CHAPTER  II 


SHOCK  CALCULATIONS 


1.  General  Solutions 

Experimental  observations  indicate  that  steady-state 
flows  exist  behind  a  shock  front  and  the  influence  of  a 
growing  boundary  layer  on  the  thermodynamic  variables  of 
the  respective  flows  is  small  within  the  observation  region. 
Therefore,  the  shock  tube  flow  may  be  treated  as  a  one  di¬ 
mensional  steady-state  problem  in  a  non-viscous  fluid — an 

12 

approximation  which  promises  the  merit  of  simplicity  . 

The  transition  from  a  state  of  low  pressure  and  temp¬ 
erature  to  one  of  high  pressure  and  temperature  is  governed 

by  the  conservation  laws  of  mass,  momentum;  and  energy  across 

13  * 
the  shock  front  .  Following  the  recent  general  formulation 

14 

outlined  by  Kim  ;  if  the  primary  and  reflected  shocks  are 
numbered  in  order  1,  2,  3>  etc.  then  for  the  n  shock  frv/ui/ 
we  have, 

(Mass)  ^On.(  v)W  r=  jOn  vK  ,  (1) 

o  2 

(Momentum)  /O,.,  vn-(  +  P*-,  =  /On  v*  +  P*  ,  (2) 

2  2 

(Energy)  hh„  +  ^v*.,  =  hh  +  |vn  ,  (3) 

where  the  subscripts  n  refer  to  variables  behind  the  n 


*  The  general  formulation  is  no  longer  than  that  for  the 
primary  shock  alone.  Formulation  is  detailed  in  order  to 
demonstrate  Kim's  new  method  completely. 
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shock  and  n-1  to  variables  ahead  of  the  shock. 

The  specific  enthalpy,  h,  is 

h  =  e  +  p^o  ,  (4) 


and  ,  v,  p  and  e  are  respectively  the  density,  the  flow 
velocity  in  the  local  coordinate  system  moving  with  the  shock 
front,  the  hydrostatic  pressure,  and  the  specific  interna?., 
energy. 

The  equation  of  state  of  the  gas  is 
p  =  RT  ,  or  p  NkT  ,  (3) 

P  r  P 


where  N,  k,  T,  R  and  f  are  the  specific  particle  number,  the 
Boltzmann  constant,  the  absolute  temperature,  and  the 
molecular  weight  of  the  gas.  Since  the  plasma  is  a  mixture 
of  different  particles  we  use  the  effective  molecular  weight 


r 


r±  h, 


i 


where  (Tj_  is  the  molecular  weight  of  the  1th  spec '.es 

and  is  its  fractional  abundance  (by  particle  number). 

For  convenience  define  Xn  and  Yn,  the  density  ratio  and 

pressure  ratio  across  the  shock  front  as 

XK  =  P^_  and  =  p* _  .  (6) 

*«-> 

Now,  (6)  in  (1)  gives  v*  =  V,  /xn  .  (7) 

(6)  and  (7)  in  (2)  gives  vh2,A.,(i-  Ij  =  Pn.(  (yh  -l)  , 


1  i 


V,  =  p,.  Xh(Yh-  1)  .  (8) 

(X*  -  1) 

i 

(8)  in  (7)  gives  (Yk  -  1)  .  (9) 

A-  K{\-  V 


(8)  and  (9)  in  (3)  re-arranges  as 


hn"  h* 


P*-/ 


j&n-i 


(Y, 


-  1)(Xk  +  1) 

“T3C 


( 10a) 


or  as,  i\  -  h*.,  =  i(pK  -  pK„  )/ 1  +  1 _ )  (10b) 

whirn  is  the  generalized  Hugoniot  relation  for  shock 
processes. 

The  total  specific  internal  energy  will  be  the  3ura  of 

the  total  translation  energy  (3  kT  per  particle),  the 

2 

rotation  and  vibration  energies  (^kT  per  particle,  per 
degree  of  freedom),  the  electronic  excitation  energy,  and  the 
energy  of  dissociation  for  molecules  and  ionization  for 
atoms. 

Therefore,  using  (5)  hn  =  eK  +  becomes 


h. 


+ 


+  fy 


N;A 

h 


+  p<j_ 
/* 


where  f*  .  f„  and  fe  are  the  rotation,  vibration  and  excitation 
degrees  of  freedom  for  the  1  particle  species.  X^anc*  °^i 
are  the  ionization  or  dissociation  potential  and  the  degree 
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of  ionization  or  dissociation  for  the  ith  species. 
For  simplification 

D.  =  2 


and  B. 


L 

4  tjjf,  ♦  r.  +  f.  \  \  • 


Therefore,  ha  =  \  (  ( B ^  +  1 )  Pn.  +  Dn 

v. 


(ID 


/On. 


Put  (li)  in  (10a),  multiply  by  A- 

'  p*-. 

to  get 


and  re-arrange 


Bh_,  +1  +  (D,w  -  D, 


(i  -  v-Hx,,  +  lj  . 


)  /O-H-l  -  (^  +  1)  Pk  /On- 1 

*  Pn-i  /%  '&0 


Let  K„  =  B*.,  +  (Dn.(  -  Dj  Ai 


(12) 


(13) 


and  solve  for  XK  to  got, 


1  + 


Bk  Y„ 


00 


K„  +  Yk 


which  is  the  generalized  Rankine-Hugoniot  relationship 
expressing  the  density  ratio  in  terms  of  the  pressure  ratio. 
If  (14)  is  substituted  in  (8)  for  X„  we  get. 


0(1  +  BkYJ  . 

b7y;  -  o-  yT) 


If  we  no'w  let 


C 

^ n  * 

Ptw 


(18) 


(16) 


(15)  can  be  re-arranged  as. 
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B.Y?  +  (1  -  Bj(l  +  rn)Yh  +  (K„rn  -  rA  -  1) 


a  quadratic  equation  for  Y*. .  Since  negative  Y*.  has  no 
physical  meaning 

Y"=  ^{[(B,-D20+O2  +  ^(1+r„  -K^f+ 


(B„  -  1)(1  -  Xj  >  . 


(17) 


If  X  and  Y  are  known,  then  from  (5a)  Tn 


JLXa.  • 
C, 


/On  =  /O.  (X^V.Xn.2  ...  X*X,) 
pK  =  P.  (Y*Yn.,  ...  YaY,  )  , 


Tn  —  To  (fc  P>i 


Ph  X?u 

p« 


R  tfc  P * 

/Ok 


NOW  V 


( VK  +  U*.,  )  ,  With  Ue  =  0,  where  VK  and 


Un-,  are  respectively  tne  shock  velocity,  and  the  gas  plasma 
velocity,  measured  in  the  laboratory  coordinate  system. 
Therefore,  to  find  the  conditions  behind  the  n^  shock 
front  it  is  necessary  to  measure  only  the  initial  gas  cco- 
dlticns  (p0  and  Tc )  and  all  n  laboratory  shock  front  vel¬ 
ocities,  Vn  .  The  gas  conditions  and  U,  may  then  be  calculated 
for  the  primary  shock  (n  =  1),  to  give  the  initial  conditions 
for  shock  front  #2  (reflected  shock)  from  which  the  con¬ 
ditions  behind  that  front  may  be  obtained.  This  procedure 
is  repeated  for  each  front  ..n  order,  until  the  nw  plasma 
conditions  are  known. 
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The  initial  and  final  conditions  of  the  gas  and  the 

degrees  of  freedom  affect  Xx  and  YK.  Prom  (14)  and  (17) 
it  cun  be  seen  that  XK  and  which  are  often  used  to 
measure  the  shock  strength,  are  ac  ually  f'unc  'ms  of  rA. 
rx  is  thus  a  good  measure  of  shock  strength. 

The  speed  of  sound,  c,  for  a  gas  is  given  by 


where  f  is  the  ratio  of  the  specific  heats  for  the  gas. 

Substitution  of  (18}  in  (16)  shows  rn  =  f  ,  where 
Mk  (=  v„.,  / c„,y )  is  the  shock  mach  velocity  relative  to 

plasma  n-1. 

Since  c  is  proportional  tc  (T  ,  a  gas  composed  of 

heavy  molecules  will  have  a  lower  speed  of  sound,  and  a 
higher  Mach  number  (stronger  shock)  for  the  same  labora¬ 
tory  velocity. 

2.  Computations 

If  we  know  the  plasma  conditions  ahead  of  the  nth 
shock  and  calculate  the  conditions  behind  the  shock  assum¬ 
ing  Bn  =  ^  ?nd  Dn  =  0  we  get  the  "ideal  gas"  con¬ 

ditions.  This  would  be  the  thermodynamic  conditions  attain- 

* 

ed  either  for  an  ideal  gas  ,  or  for  the  gas  immediately 

*  In  this  paper,  the  term  ideal  gas  denotes  a  gas  with  no  in¬ 
ternal  degrees  of  freedom.  Van  der  kaals  forces  may  always  be 
neglected  for  gas  conditions  encountered  in  the  shock  tube. 
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oehlnd  the  shock  front  after  relaxation  of  only  the  trans¬ 
lational  energy  states.  Other  energy  states  remain  unchang¬ 
ed  and  would  temporarily  exhibit  population  distributions 
corresponding  to  the  original  temperature.  It  is  from 
these  gas  conditions  that  relaxation  into  the  other  states 
actually  occurs.  For  this  special  case  calculations  become 
very  simple,  especially  for  the  primary  shock,  and  com¬ 
putations  can  readily  be  done  by  hand.  Full  computations 
for  the  final  equilibrium  state  of  the  non-ideal  gas  will 
require  the  use  of  a  computer. 

Relaxation  studies  are  usually  limited  to  the  primary 
shock  region  where  the  relaxation  times  are  longest  and 
where  computations  rre  somewhat  simpler.  For  the  n  =  1 
shock,  we  will  have  B0S  4  and  D  s  0  for  all  reasonable 
initial  temperatures.  The  use  of  monatomic  gases  simpli¬ 
fies  both  the  actual  processes  involved  in  relaxation 
(only  electronic  excitation  and  ionization  exist),  and  the 
computation  of  final  conditions,  since  f*  =  f*  s  0  ,  and  D, 
involves  only  ionization. 

The  calculations  for  values  of  X,  and  Y,  (hence  p,  ,p, , 

Ti )  involve  the  summation  of  each  electronic  partition 
function  (Z;)  to  some  "cut-off"  energy,  and  the  degrees  of 
ionization  (  )f  but  the  cut-off  energies  and  the  Zx'  , 
themselves  depend  on  -  he  unknown  final  gas  conditions 
(p,  ,  T,  ).  Final  conditions  may  be  assumed  In  order  to 
evaluate  the  and  ZA  for  use  in  approximating  the  final 
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conditions.  These  better  values  are  then  used  to  improve  the 
°<i  and  Z i  which  will,  in  turn,  result  in  even  better  final 
values.  This  iteration  process  may  be  repeated  until  the  in¬ 
put  and  output  values  agree  to  within  any  desired  limit. 

6  14 

Kim  claims  J  that  this  simple  iteration  process 
converges  very  rapidly,  but  in  fact  divergence  may  actually 
result  unless  the  original  choices  of  values  were  reasonably 
close  to  the  correct  values.  Rapid  convergence,  even  when 
the  original  choice  is  quite  poor,  may  be  easily  obtained 
if  some  form  of  interpolation  based  upon  previous  input- 
output  values  is  incorporated  after  each  iteration  cycle, 
to  obtain  even  better  values  for  use  in  the  next  cycle. 


For  a  plasma  of  electrons,  atoms  and  ions  we  may  write 
the  Saha  equations  for  the  first  two  ionizations  of  the 


atoms  as 


N+  Nc 

N, 


Z+  Z*  e 


"tyr 


(19) 


and  N++  N« 
N* 


Zi-t  Ze 
Zy 


'W 


(20) 


Hi  >  Ne ,  N+  and  are  the  number  densities  per  unit  volume 
for  atoms,  electrons,  singly  and  doubly  ionized  atoms,  with 
ZA  ,  Ze,  Z*  and  Z+t  as  the  corresponding  partition  functions, 
and  are  the  effective  Is ^  and  2nd  ionization  poten¬ 


tials  of  the  atom. 
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We  now  define  the  1st  and  2nd  degrees  cf  Ionization  as 


Prom  which  if  Ne  =  N„  +  +  N,+  we  have, 

N„  =  Ne(l  -<*,),  Nv  =  Nc  ( 1  -  o<A )  cxCt  , 

Nf+  =  N0(^,p<J,  and  Ntf  =  <*,  N0  ( 1  +  «<, ) . 

(22) 


Use  of*  p  =  NkT  and  N  =  N„  (1  +  ^  +  oC,^ ) 
with  (22)  causes  (19)  and  (20)  to  simplify  as 


Or,  o^(l  -  +  f,  +  f,  «<*  )  -  -  1*,  =0  (25) 

and  o^o(j(  1  +  fjj )  +  ( **»,  +  fj  )  “  ( 1  +  c^, )  fa  =  0  . 

(26) 

If  f,  and  f*  are  known  (or  estimated)  oCt  and 

may  be  found  quickly  by  a  simple  iteration  process  for  the 
positive  roots  of  these  quadratic  equations. 

*  All  N  now  refer  to  the  particle  number  density  per  cc. 
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If  gas  mixtures  are  used  the  same  method  is  employed  but 
further  Saha  equations*  and  a  revised  N  are  necessary  to  in¬ 
clude  the  atoms  of  the  second  gas.  The  treatment  of  mixtures 
with  polyatomic  gases  has  been  given  by  Wilkerson^  and  was 
used  only  briefly  for  this  research  in  connection  with  the 
impurity  studies. 


For  any  species  in  the  plasma,  the  partition  function  is 
given  by  Z  =  ZrZRZvZ/r  with  the  index  T  standing  for 
translation,  etc.  The  atoms  and  ions  are  identical  except 
for  charge  and  very  small  differences  in  mass,  so  only 
their  Ze  may  differ.  For  electrons 

Zr  =  /  2T  mkT  fa  ,  Z„  =  2  (2  spin  states) 

1  h2  / 

and  Zv  =  1  . 


Using  these  values  in  (23)  and  (24)  for  f,  and  f 


2T  m 
h2 


,  -Mr 

)  (kT)  Ze+  e 
/  P  Ze,, 


a  > 


(27) 


and 

h  = 


Mf  Z£«_  e'^7-  . 

P  Zet. 


(28) 
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It  has  been  seen  experimentally  that  for  high  pressures 
and  charge  densities  atomic  spectra  remain  normal  up  to  a 
certain  quantum  number  after  which  the  lines,  merge  Into  n 
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continuum  with  the  spectral  series  limits  depressed  as 
though  a  lower  ionization  potential  existed. 


Thus 


(29) 


where  gnJ  E„,  and  are  the  statistical  weights,  the  ex¬ 
citation  energies  and  the  effective  (lowered)  ionization 
energy  for  the  atom  or  ion  considered.  All  gn,  En  and 

the  free  space  ionization  potentials  for  computation  were 

17 

obtained  from  the  tables  by  Charlotte  Moore  1 . 
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Inglis  and  Teller  took  stark  broadening  as  the  impor¬ 
tant  mechanism  for  level  spreading  and  found  the  cut-off 
quantum  number  (or  )  by  equating  the  energy  difference 
of  neighboring  levels  to  the  spread  of  the  corresponding 
levels.  Better  results  take  the  charge  screening  of  the 
plasma  into  account  as  well^'^.  The  Debye-Huckel  formula 
for  the  lowering  of  the  ionization  potential  (eq.  30)  is 
obtained  this  way  but  is  valid  only  when  the  number  of 

electrons  per  Debye  sphere  is  less  than  one.  For  higher 
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density  plasmas  the  result  of  Ecker  and  Kroll  should  be 
used. 

For  the  present  work  the  Debye-Huckel  cut-off  is 
satisfactory.  Thus, 


=  e2  (z  +  1)  (30) 

V/j 

__j<t _  ,  (3D 

(Ne  +  H  z2N;  )  j 

z 


with 
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where  z  =  0,1,2  etc  for  neutral  atoms,  1st  ions,  2nd 

ions,  etc. 

The  actual  ^  used  here  is  relatively  unimportant, 
but  some  consistent  cut-off  procedure  is  necessary  other¬ 
wise  all  are  infinite.  Tests  show  that  even  for  rela¬ 
tively  strong  shocks  tripling  A\z  causes  less  than  1 % 
variation  in  the  final  thermodynamic  variables  calculated. 

If  the  final  p,  T  are  known,  or  assumed,  the  A  may 

be  obtained  from  (31)  for  each  species.  These  values  allow 
the  effective  to  be  found  for  use  in  (27),  (28)  and 

(29).  f,  and  f2  can  then  be  found  for  use  in  (25)  and  (26) 
where  and  are  found  by  iteration.  The  shock  equations 
(1*0  and  (17)  may  then  be  solved  to  find  better  values  for 
the  final  p,^o  and  T  which  in  turn  are  used  (with  inter¬ 
polation)  to  begin  the  next  cycle  of  iteration. 


CHAPTER  III 


CONSTRUCTION  OP  THE  SHOCK  TUBE 


1.  Introduction 

A  heated  shock  tube  has  been  in  full  operation  at  the 

University  of  Michigan  since  1967.  Designed  and  built  by 

Y.W.  Kim,  its  initial  application  was  the  study  of  high 

electron  density  mercury  plasmas.  Preliminary  observations 
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and  construction  were  reported  earlier  3 

During  early  work  by  Kim,  some  minor  design  difficulties 
were  noted,  and  this  chapter  describes  these  problems  and 
recent  engineering  Improvements  to  correct  them.  Some 
suggestions  for  future  heated  tubes  are  included.  A  more 
complete  description,  which  includes  details  for  a  new 
flush  valve,  and  the  construction  for  large  (7")  observation 
windows,  may  be  found  in  reference  23.  Fig.  2  shows  the  en¬ 
tire  shock  tube  schematically. 

The  main  advantages  of  a  heated  tube  are  the  high  purity 
levels  of  the  shocked  gas,  resulting  from  the  continual 
baking,  and  the  ability  to  work  with  substances  evaporable 
only  at  moderate  temperatures. 

2.  The  Expansion  Tube.. 

The  main  shock  tube  with  connected  fittings  and  hardware, 
was  constructed  of  stainless  steel  (Type  30*0.  The  expansion 
tube,  Including  the  test  section,  is  107  long  and  consists 
of  five  sections,  made  from  seamless  rectangular  tubing 


with  3/8"  wall,  inner  limensions  1.510"  x  2.560",  and 

inner  corner  radius  0.300". 

In  order  to  reduce  boundary  layer  effects,  the  inner 

surface  of  thr  r-ntire  expansion  tube  was  highly  polished. 

Tube  sectioning  was  considered  necessary  because  the 

rectangular  cross-section  is  extremely  difficult  to  polish 

in  long  sections.  The  sections  were  put  together  by 

flanges  7/8"  thick,  and  6"  in  diameter.  Between  each  pair 

of  adjoining  sections,  goes  a  6"  diameter,  1-3/4"  thick 

disc,  containing  two  viewing  or  pro be -mounting  ports  fitted 

for  2-3/4"  Conflat  flanges.  Details  may  be  seen  in  Pigs. 

3  and  4.  The  matching  faces  of  the  disc  and  adjoining 

flange  have  gasket  grooves  similar  to  those  of  the  6" 
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Varian  Conflat  flanges 

OFHC  copper  gaskets,  0.080"  thick,  are  sealed  between 
the  matching  conflat  grooves  using  ten  3/8" -?4  x  4"  bolts, 
with  the  flanges  tightened  so  they  come  into  direct  face- 
to-face  contact.  Aligning  is  done  by  means  of  preset 
dowel  pins,  and  the  sections  are  not  interchangable.  All 
nuts  and  bolts  are  lubricated  with  Sauereisen  graphite 
paste  (No.  51),  to  prevent  them  from  fusing  together  at  the 
high  temperatures. 

Kim  found  that  after  a  period  of  use  any  change  in  tube 

temperature  caused  a  small  vacuum  leak  at  one  of  the  many 

flanges.  This  could  often  be  corrected  by  further  bolt 
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tightening  on  the  correct  flange  '  .  Normally,  neither  of 
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these  should  happen,  according  to  Wheeler  and  Coulder  , 


Figure  4.  Partially  Assembled  Heated  Shock  Tube 


who  researched  the  Conflat  type  flange  extensively.  However, 
their  paper  did  suggest  several  necessary  modifications. 

First,  the  Conflat  seal  is  due  to  the  interaction  of  the 
outer  circumference  and  the  knife-edge  only,  and  the  inner 
part  of  the  gasket  should  be  free  from  any  compression  or 
restriction.  Second,  the  nickel  plating,  which  was  used 
to  protect  the  gasket  from  the  mercury  vapor,  may  chip  or 
crack  at  the  knife-edge,  thus  causing  leaks.  Third,  the 
knife-edge  should  cut  into  the  gasket  0,015"  on  each  side 
for  best  results.  When  the  tube  was  dismantled  it  was  found 
that  'he  knife-edges  had  been  badly  blunted  by  the  harder 
plated  gaskets.  These  three  faults  were  corrected  by  re¬ 
machining  all  flanges,  and  eliminating  gasket  plating. 
Overhauls  will  now  be  necessary  about  every  two  years  since 
the  unplated  gaskets  are  more  rapidly  attacked  by  the  heated 
mercury  vapor. 

A  fourth  problem,  which  probably  aggravated  the  effect 

of  the  other  three,  and  vice  versa,  was  the  asymmetric  cool- 

pa 

ing  of  the  flanges.  This  has  definitely  been  noted  to 
cause  leaks.  The  expansion  tube  was  originally  supported  by 
a  pair  of  short,  large  cross-section  legs  at  each  flange 
(see  Fig.  4).  To  correct  this  heat  loss  problem,  the  legs 
were  removed,  and  the  tube  suspended  from  the  outer  casing 
by  four  pairs  of  6"  long,  1/4"  ctainless  steel  eye  bolts, 
fastened  to  the  smaller  observation  port  flanges  on  the  discs. 
In  designing  a  new  tube  this  suspension  would  be  best,  but 
with  the  eye  bolts  attached  t)  special  brackets  welded  to 
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the  tube  between  flanges.  With  such  a  suspension  an  extra 
1/4"  diameter,  6"  bolt,  between  the  end  of  the  test  section 
and  the  outer  casing, was  necessary  to  prevent  tube  recoil. 

Also,  since  each  disc  requires  6  seals  which  must  be 
maintained,  (2  large  and  2  small  Conflat,  and  2  windows) 
it  would  be  wortn  the  effort  to  polish  the  tube  in  longer 
pieces.  It  is  also  worth  considering  the  use  of  a  circular 
shock  tube,  which  can  be  cheaply  and  easily  honed  as  one 
piece,  but  flow  studies  would  then  be  difficult  to  inter¬ 
pret.  To  cut  heat  losses,  machining,  and  leaks,  the  discs 
should  be  eliminated  completely  and  the  observation  ports 
mounted  directly  on  the  tube  wall  using  a  structure  similar 
to  that  described  for  the  newest  test  section. 

3.  The  Compression  Tube 

The  compression  section  is  constructed  as  a  single 
piece  48"  long,  of  the  same  rectangular  tubing.  Flanges  1" 
thick  are  welded  on  each  end  between  the  tubing  and  a 
surrounding  steel  tube,  6"  In  diameter  and  1/4"  thick,  and 
the  3pace  between  filled  with  Sauereisen,  low  expansion, 
high  temperature  cement  (No.  7). 

The  compression  tube  is  connected  to  the  expansion  tube 

by  six  sets  of  5/8" -10  x  4",  heat-treated  nuts  and  bolts. 

A  round,  soft  metal  diaphragm  is  squeezed  between  the  faces 

of  matching  flanges.  Two  concentric  knife-edges  on  each 

flange  complete  the  vacuum  -  high  pressure  packing.  The 

-5 

total  seal  proved  good  from  5  x  10  mm.  Hg.  to  3000  psi. 
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For  replacing  the  diaphragm  after  each  shot,  the  com¬ 
pression  tube  rolls  on  two  pairs  of  Fafner  ball-bearing 
casters  mounted  on  a  slotted  angle  rail  structure  in  a 
section  of  the  outer  casing. 

4.  Heating 

Although  it  performed  adequately,  Kim's  original  cylin¬ 
drical  oven^  occasionally  developed  electrical  shorts  causing 
the  wire  winding  to  burn  out.  Also,  its  weight  and  bulk  were 
a  severe  handicap  during  maintenance  and  overhaul.  These 
faults  were  overcome  with  a  new  heating  structure  which  also 
allows  much  more  rapid  heating  and  cooling  of  the  shock  tube, 
and  requires  only  one-tenth  the  power  consumption.  (100  watts 
per  foot  at  500 °C). 

The  simple  shock  tube  sections  are  wrapped  with  4"  wide 
asbestos  tape,  then  #16  Chromel  "C,:  heater  wire  is  wound  over 
this,  and  the  structure  covered  by  Sauereisen  electric  heater 
cement  (No.  63).  The  test  section  and  the  mercury  reservoir 
area,  being  complicated  shapes,  are  heated  by  means  of 

"Hot-Foil"  heating  tapes.  The  current  to  each  section  is 
controlled  individually  by  variac  to  give  uniform  tube  tem¬ 
perature,  checked  by  thermocouples.  Since  the  temperature 
of  the  tube,  once  established,  was  very  stable  no  special 
automatic  temperature  control  was  needed. 

The  Chromel  wire  is  far  heavier  than  necessary  for  the 
steady  state  currents  used  even  at  5°0°C,  thus  allowing 
current  Increase  for  rapid  heating  during  temperature  changes 


without  wire  burn-out.  The  cement  and  asbestos  serve  as  elec¬ 
trical  insulation  from  the  tube.  The  cement  also  protects  the 
asbestos  from  air  since  its  binding  material  deteriorates  at 
very  high  temperature,  and  in  the  future,  if  necessary,  will 
hold  the  decomposed  asbestos  and  the  wire  in  place. 

Two  layers  of  1^"  thick  fibi^glas  pipe  covering  (7"  and 
11"  inner  diameters)  were  fitted  around  the  shock  tube  inside 
the  outer  casing.  Joints  were  staggered  to  prevent  air  con¬ 
vection  to  the  outer  casing  as  much  as  possible.  The  highest 
temperature  usable  is  determined  by  the  deterioration  of  this 
insulation  and  should  be  in  excess  of  500° c.  More  insulation 
is  pres  -i  on  the  outside  of  the  casing.  The  entire  area  in¬ 
side  the  casing  is  not  filled  by  insulation  because  the  uni¬ 
formity  of  the  tube  temperature  depends  on  free  air  convect¬ 
ion  near  the  tube. 

Iron-constantan  chermocouples  were  attached  at  6  points 
on  the  tube.  The  test  section,  the  reservoirs,  the  4  disc 
pieces  and  each  end  of  the  compression  tube  are  monitored 
points.  The  thermocouples  are  resistance  matched,  and  using 
a  multipole  switch  are  checked  individually  by  connecting 
the  desired  thermocouple  to  a  0  -  260° C  API  Pyrometer 
(Model  603-1).  Temperature  can  be  controlled  to  within  |°C. 

5.  The  Outer  Casing 

The  outer  casing  is  also  divided  into  sections,  with 
lengths  roughly  matching  those  of  the  shock  tube,  and  made 
of  half-hard  temper  brass  tubing  with  a  14"  inner  diameter 
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and  3/16”  thickness.  The  sections  are  joined  by  firebox 
quality  steel  flanges.  Wrapped  around  each  tubing  is  1^" 
flberglas  insulation  covered  with  4"  wide  asbestos  tape 
making  the  overall  diameter  almost  18".  Each  section  has 
two  side  openings  of  5-1/4"  diameter,  allowing  access  to  the 
shock  tube  at  the  viewing  port  positions.  One  12"  diameter 
opening  at  the  middle  section  accommodates  all  the  vacuum 
components  to  the  main  tube  and  mechanical  feed-throughs 
to  the  mercury  reservoirs. 

The  failure  of  an  observation  window  or  of  other  pack¬ 
ings  of  the  shock  tube  may  cause  exposure  of  laboratory 
personnel  to  the  extremely  hazardous  mercury  vapor.  The 
present  safety  system  depends  upon  the  casing  being  moderate 
ly  air-tight.  The  inside  of  the  casing  is  connected,  at  the 
test  section  end, to  a  long  diameter  pipe  opening  directly 
outdoors.  This  prohibits  any  gas  build-up  inside  the  casing 
which  would  cause  leaks  into  the  laboratory.  When  working 
with  poisonous  vapors  the  outer  casing  is  never  opened 
after  firing  a  shock  until  the  high  pressure  is  released, 
and  an  attempt  is  made  to  evacuate  the  shock  tube.  If  a 
good  vacuum  is  not  obtained,  then  the  casing  is  flushed  with 
compressed  air  to  drive  any  remaining  vapors  from  the  cas¬ 
ing. 

With  the  new  improvements  the  brass  casing  now  serves 
no  useful  purpose,  in  fact  it  is  a  handicap  during  any 
repair  or  overhaul.  A  new  "casing"  is  suggested  for  future 
tubes,  which  is  cheaply  and  easily  constructed.  The  shock 
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tube  is  suspended  by  longer  eye  bolts  to  a  metal  frame,  but 
with  extra  bars  fastened  near  the  diaphragm  flanges  to  pre¬ 
vent  recoil  motion.  Three  layers  of  fiberglas  pipe  covering 
surround  the  tube  (7",  10"  13") ,  with  joins  staggered,  and 
wrapped  with  4"  wide  asbestos  tape.  This  "casing"  is  lined 
with  aluminum  foil  as  a  vapor  barrier,  either  inside  or  be¬ 
tween  the  layers  of  pipe  covering,  and  suspended  from  the 
metal  frame.  Electrical  and  instrumental  access  may  be  ob¬ 
tained  directly  through  the  insulation  without  causing 
appreciable  leaks. 

Optical  clearance  for  a  small  observation  port  is  easily 
obtained  through  a  thin -walled  steel  tube  which  is  in  con¬ 
tact  with  the  port  at  one  end  and  passes  through  the  layers 
of  insulation  to  the  outside.  This  viewing  system  is  actual¬ 
ly  used  with  the  new  test  section.  A  light  weight  window 
frame  can  easily  be  fitted  into  the  insulation  to  clear  an 
area  for  larger  test  section  windows.  The  new  heating, 
thermocouple,  and  safety  systems  are  still  applicable.  The 
insulation  casing  for  the  compression  section  may  be  fitted 
tightly  to  the  tube  arid  the  entire  section  suspended  from 
a  carriage  which  rolls  back  and  forth  on  the  metal  frame. 

6.  The  New  Test  Section 

A13  window  seals  of  the  heated  shock  tube  must  be  highly 
vacuum  and  pressure  tight.  The  heating  requirements  and  the 
presence  of  mercury  specify  the  use  of  stainless  steel  0- 
rings  and  soft  metal  gaskets.  These  are  adequate  for  metal- 
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to-metal  seal3,  but  for  optical  windows  metal  O-rlngs  tend 
to  chip  the  surface,  and  fusing  metal  into  optically  flat 
glass  is  not  easy.  Instead,  direct  contact  between  the 
glass  and  metal  O-rlng  Is  avoided  by  using  a  thin  layer  of 
Kapton  polymlde  film,  called  H-film.  The  H-film  is  avail¬ 
able  in  0.001"  to  0.010"  thicknesses,  and  has  proved  highly 
satisfactory. 

H-film  is  extremely  strong  and  retains  its  strength  over 
a  temperature  range  far  greater  than  used  here.  Because  H- 
film  does  not  "flow"  extensively  under  pressure  the  best 
seal  is  made  by  pressing  an  area  only  about  1/16"  wide, 
about  the  area  under  a  stressed  metal  0-rlng.  Broader  areas 
require  increased  bolt  loads,  and  often  cause  window  crack¬ 
ing  due  to  increased  stresses. 

Vacuum  seals  for  large  windows  are  extremely  difficult 
to  achieve  and  maintain  at  high  temperatures.  After  moder¬ 
ate  success  with  attempts  to  improve  Kim's  original  large 
(7")  window  test  section,  it  was  decided  that  the  remaining 
faults23  could  only  be  corrected  by  building  a  costly  new 
7"  test  section.  Instead,  a  very  simple  structure  was  tried 
which  proved  to  be  excellent,  and  solved  all  the  window 
problems. 

The  new  test  section,  shown  in  Fig.  5,  consists  of  a 
regular  piece  of  the  shock  tube  with  3  observation  ports, 
or  windows,  2.750"  apart  (at  20° (;)  along  each  side,  and  the 
final  window  1.75"  from  the  end-wall.  Windows  on  opposite 
aides  are  exactly  aligned  so  that  beams  of  light  may  be 


Figure  5.  Cutaway  Top  View  of  New  Test  Section. 
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passed  through  the  tube  perpendicular  to  the  tube  axis. 
Details  of  the  window  structure  are  shown  in  Pig.  6. 

A  vertical  arrangement  of  3  holes  is  used  in  each  win¬ 
dow  to  allow  more  light  to  be  collected.  Holes  are  used, 
instead  of  a  narrow  vertical  slit,  for  minimum  gas  flow 
disruption.  The  windows  are  set  into  the  wall,  and  the  hole 
size  stepped,  to  keep  reflections  from  hole  walls  to  a  min¬ 
imum.  The  raised  sealing  ring  allows  all  pressure  to  be 
limited  to  a  l/l6"  wide  area  on  the  H-film,  and  removes 
pressures  from  the  edge  of  the  glass.  O-rings  are  not  used 
because  they  are  difficult  to  change  after  use,  and  the  0- 
ring  groove  would  weaken  the  tube  wall. 

The  glass  windows,  3/8"  thick  and  1"  in  diameter,  may 
be  replaced  by  quartz  or  other  material,  if  desired. 

7.  The  Mercury  System 

A  special  gauge,  and  a  mercury  reservoir  are  installed 
within  the  oven,  close  to  the  shock  tube.  This  system  must 
be  kept  at  the  same  temperature  as  the  tube,  or  higher, 
since  the  maximum  vapor  pressure  usable  will  be  that  of  the 
coolest  part  of  the  entire  system.  The  gauge  and  reservoir 
are  separated  from  one  another,  and  connected  to  the  tube 
using  two  bakable  valves  made  with  stainless  steel  bellows 
and  0-rings.  The  arrangement  is  shown  in  Pig.  7. 

The  common  difficulties  encountered  in  measuring  pres¬ 
sures  of  heated  vapors  are  the  condensation  at  cold  spots 
caused  by  connecting  a  gauge  at  roo a  temperature,  and  the 
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Tube  Wall  ^  0.060" 


A.  Top  Cutaway  View  (Double  Actual  Size) 


Figure  6.  New  Teat  Section  Window  Structure 
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failure  of  many  molecular  gauges  when  excessively  heated. 

As  a  result,  Y. W.  Kim  designed  a  gauge,  shown  in  Pig.  7 
which  can  operate  at  the  tube  temperature. 

The  space  enclosed  by  a  stainless  steel  sylphon  bellows 
(spring  rate  17.9  lb/in.  at  20°C. )  is  filled  with  the  gas  to 
be  measured  and  the  bellows  are  housed  in  a  small  vacuum 
cylinder  so  that  the  pressure  outside  the  bellows  may  be 
independently  regulated.  This  method  eliminates  calibration 
which  would  be  necessary  if  the  bellows  expansion  were  used 
as  the  pressure  scale --the  spring  rate  varies  with  temper¬ 
ature.  The  effect  of  cooling  and  of  the  temperature  gradient 
on  the  pressure  distribution  in  the  outer  system  are  neglig¬ 
ible. 

An  insulated  electrode  is  sealed  into  the  outer  cylinder 
and  extends  far  enough  to  make  contact  with  the  top  of  the 
bellows  when  both  sides  are  fully  evacuated.  A  microammeter 
is  used  to  measure  the  contact  current.  Before  any  gas 
measurement,  the  air  pressure  outside  of  the  bellows  is  raised 
until  the  electrode  contact  is  nearly  broken,  at  say  3°  pAmp. 
The  air  pressure  is  then  measured  by  an  oil  manometer  and 
used  as  a  zero  point.  The  air  pressure  is  then  raised  by  the 
desired  mercury  gas  pressure  (meter  reads  zero  now),  and  the 
mercury  gas  is  then  filled  into  the  bellows  (and  shock  tube) 
from  the  liquid  reservoir  until  the  gauge  contact  is  re¬ 
stored.  The  outside  air  pressure  is  then  carefully  adjusted 
until  the  meter  again  reads  30  p  amp.  The  final  manometer 
leading  corrected  by  the  zero  point  vaiue  gives  the  mercury 
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Figure  7.  Bellows  Gauge  and  Mercury  Reservoir 
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gas  pressure. 

This  filling  and  measuring  procedure  allows  accurate 
determination  of  the  pressure  (but  not  an  exact  pre-determin- 
ed  value)  while  keeping  impurity  contamination  to  a  minimum. 
The  liquid  mercury  in  the  reservoir  was  thoroughly  outgased 
before  use  to  guarantee  purity,  but  tube  leaks  and  outgasing 
are  always  present.  Experimental  shocks  were  never  fired 
when  the  total  tube  leak  (plus  outgasing)  was  greater  than 
0.2  ju/minute,  which  gives  a  maximum  of  0.6  p.  impurity  in  the 
3  minutes  necessary  to  close  and  fill  the  tube  before  firing 
the  shock.  For  the  lowest  values  of  Po  used  (10  Torr)  this 
0.6  p  represents  an  impurity  of  60  parts  per  million. 

Maintenance  of  this  bellows  gauge  is  simple  and  the 
stainless  steel  sylphon  bellows  have  very  long  life.  The 
present  bellows  has  been  in  use  since  construction  3  years 
ago.  However,  the  electrical  contact  within  the  gauge 
sometimes  becomes  bad,  and  requires  opening  the  gauge  to 
polish  the  surfaces  involved. 

8.  Gas  Handling  Systems 

A  Duo-Seal  pump  (Type  1  402-B)  backs  a  CVC  metal  mercury 
diffusion  pump  (Type  MHG-180)  to  evacuate  the  tube  and  the 
gas-vacuum  complex.  A  CVC  glass  trap  between  the  pumps 
guards  against  mercury  contamination.  The  diffusion  pump  is 
connected  to  a  manifold  through  a  stainless  steel  water 
cooled  trap,  and  the  bakable  flush  valve  connects  the  mani¬ 
fold  to  the  shock  tube  and  closes  t.ie  tube  off  completely 
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from  the  vacuum  system  when  a  shot  is  fired.  Also  connected 
to  the  manifold  are  an  Autovac  gauge  (Type  3294-B),  and  an 
oil  manometer,  which  are  employed  for  monitoring  the  various 
pressures. 

Within  one  minute  the  expansion  section  is  ev  ^cuated  to 
1  p  Hg  and  ten  minutes  pumping  brings  the  tube  pressure  to 
0.001  p  Hg,  with  a  leak  rate  of  about  0.1  p/ minute,  brief 
exposure  of  the  tube  to  ai**  has  no  significant  effect  on  the 
pumping  rate,  however,  an  hour  is  required  to  bring  the  tube 
back  to  uniform  temperature  after  the  outer  casing  is  opened 
for  a  diaphragm  change. 

Typical  gases  used  to  drive  shocks  are  hydrogen,  helium 
and  nitrogen.  These  are  commercially  available  in  bottles 
at  pressures  up  to  3500  psi.  These  bottles  are  connected, 
via  stainless  stt?el  flexible  hoses,  to  the  compression 
section,  and  a  manifold,  which  allows  the  various  gases  to 
be  used  separately,  or  as  mixtures. 


CHAPTER  IV 


OBSERVATION  TECHNIQUES 

1.  The  Schlieren  System 

The  arrangement  for  the  schlieren  system  at  each  window 
of  the  test  section  is  shown  schematically  in  Figure  8.  A 
3^  watt  tungsten  filament  lamp  (G.E.  type  BXJ)  illuminates 
a  0.010"  wide  slit,  SI  as  the  source  for  3  schlieren  systems. 
Lens  LI  with  slit  SI  at  the  focus  results  in  a  parallel 
beam  of  light.  Two  sets  of  adjustable  mirrors  separate  this 
beam  into  3  parts  and  align  each  of  these  parallel  beams 
through  a  different  test  section  window.  The  emerging  beam 
passes  through  lens  L2  and  forms  an  image  of  SI  at  the  knife- 
edge  position,  the  focus  of  L2.  The  moveable  knife-edge  is 
adjusted  over  the  image  of  Si  so  that  very  little  light 
normally  reaches  S2,  a  0.006"  wide  slit  in  front  of  the 
photomultiplier  tube. 

A  density  gradient  anywhere  within  the  parallel  section 
of  the  beam  will  deflect  part  of  the  beam  and  change  the 
amount  of  light  reacting  S2.  Density  gradients  in  the  shock 
tube  may  thus  be  detected,  while  a  new  constant  density  in 
the  tube  will  cause  no  change  in  the  schlieren  light  signal. 
Best  results  are  obtained  if  the  test  section  is  imaged  at 
S2  with  unit  magnification,  (distances  =  2f).  This  is  a 
compromise  between  spatial  resolution  and  light  intensity 
at  the  photomultiplier.  The  knife-edge  is  positioned  so  that 
the  Increase  In  density  ;>a  the  nri  ary  shoot  passes  will 
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cause  an  increase  in  the  light  reaching  S2.  Note  that  small 
adjustments  of  the  knife-edge,  which  is  on  the  axis  of  the 
optical  system,  greatly  affect  the  amount  of  luminosity 
observed  at  S2. 

Ab  the  shock  front  passes  the  window,  the  gradient  is 

so  steep  and  short  that  light  passes  very  strongly  for  only 

a  fraction  of  a  microsecond.  A  fairly  insensitive  schlieren 

system  is  used  in  order  that  only  the  large  density  gradient 

at  the  shock  front  can  he  seen.  A  measure  of  the  sensitivity 

2  6 

of  a  schlieren  system  is  given  by  A I  «<  f*.  , 

- r  - gj- 

where  f2  =  focal  length  of  the  lens  L2,  d2  =  the 
width  of  the  SI  image  at  the  knife-edge,  and  «<  =  the 
angle  which  the  light  beam  is  refracted.  This  equation  is 
valid  only  for  gradients  which  do  not  move  the  Si  image  en¬ 
tirely  off  the  knife-edge. 

The  shock  front  velocity  is  determined  from  the  time 
between  the  schlieren  pulses  from  different  windows  and  the 
inter-window  distances. 

2.  Luminosity  Observation 

Because  the  test  section  is  focussed  on  S2  the  photo¬ 
multiplier  will  automatically  record  the  plasma  luminosity 
as  well  as  the  light  from  the  schlieren  system.  The  photo¬ 
multiplier  does  not  record  all  the  plasma  radiation  because 
the  glass  windows  and  tube  envelope  absorb  ultraviolet  with 
wavelength  shorter  than  3200A° ,  and  the  931 A  lube  la  not 
sensitive  to  wavelengths  longer  than  6b00A° . 
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The  circuits  for  the  photomultiplier  tubes  are  shown 
schematically  in  Pig.  9.  The  circuit  for  tube  #1  (re  shock 
passage)  is  designed  mainly  for  observing  the  plasma 
luminosity  without  distortion.  The  #2  and  #3  tubes  are  used 
only  for  shock  speed  timing  and  their  outputs  are  limited 
to  1  volt  for  pulses  longer  than  5  Msec,  duration.  The 
output  of  tube  #1  is  free  from  distortion  up  to  30  volts  for 
signals  from  0.1  Msec,  to  0.01  sec.  Output  limiting  and 
distortion  are  controlled  by  the  choice  of  resistances 
and  capacitances,  and  the  values  used  are  shown  with  Figure  8. 

A  Tektronix  Model  551  dual  beam  oscilloscope  was  used  to 
record  results.  The  output  of  photomultiplier  #1  (schlieren 
plus  luminosity)  was  displayed  on  the  lower  beam,  while  the 
outputs  of  #2  and  #3  were  added  by  the  oscilloscope  pre¬ 
amplifier  (type  CA)  and  shown  Inverted  on  the  upper  beam. 

Sweep  was  triggered  (Internally)  by  the  schlieren  pulse  on 

#1. 

The  sweep  speed  was  chosen  to  give  accuracy  in  deter¬ 
mining  the  shock  velocity.  The  gain  on  the  upper  beam 
pattern  must  be  kept  small  to  ensure  that  the  pattern  is 
not  lost  off-scale.  This  is  the  main  reason  for  limiting 
#2  and  #3  tube  outputs  to  1  volt  each.  Observations  are 
virtually  useless  if  the  shock  speed  is  not  known  (and 
hence  the  gas  conditions). 

For  accurate  determination  of  the  relaxation  time  t  , 
it  is  desirable  to  have  the  #1  pattern  large,  so  that 
changes  in  slope  are  readily  apparent,  and  to  have  ^  occupy 
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a  large  fraction  of  the  10  cm  trace.  The  output  of  tube  #1 
was  also  observed  using  a  Tektronix  Type  535A  oscilloscope 
with  the  sweep  and  gain  set  to  give  a  good  relaxation 
profile  trace. 

The  oscilloscope  patterns  are  recorded  on  type  47 
Polaroid  film,  and  a  "zero"  trace  recorded  as  well  for  ref¬ 
erence.  Contact  prints  of  typical  oscilloscope  patterns 
for  mercury  are  shown  below. 


e 


,  * 

M-4l 

Oscilloscope  535A 
Sweep  =  5  us/cm 
p,  =  29.1  Torr 

T„  =  514.2  °K 


I  ' 

t 


M-10 

Oscilloscope  551 
Sweep  =  20  us/cm 
p„  =  16.1  Torr 

T0  =  513.2  °K 


*  M  denotes  shocks  in  mercury,  A  shocks  in  argon. 
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3.  Limits  on  Observations 

The  two  pictures  below  were  obtained  for  the  same  weak 
shock  in  argon.  Neutral  density  filters  are  used  to  cut  the 
luminosity  suitably.  No  schlieren  light  or  knife-edge  is 
present  so  triggering  occurred  for  the  first  oscilloscope 
when  the  luminosity  became  large  enough,  and  the  second  was 
triggered  simultaneously  by  the  first. 


10$  of  light  from 
window  #1 

Sweep  =  20  us/cm 

Gain  =  0.5  v/cm 


0.01$  of  light  at 
window  #2 

Sweep  =  20  us/cm 

Gain  =0.5  v/cm 


These  pictures  show  examples  of  some  problems  encountered 
with  the  present  observation  technique.  In  the  first  picture, 
the  section  A  to  B  is  similar  to  the  luminosity  curve  of 
M-41.  7^  Is  very  long  and  relaxation  appears  to  be  just 

complete  when  the  interface  arrives,  at  70  Msec.  (b). 
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Longer  could  not  be  observed.  Instead  of  the  luminosity 
dropping  abruptly  when  the  interface  arrives,  a  spike  appears 
with  the  interface.  (B-*C).  The  spike  is  probably  due  to 
extra  luminosity  caused  by  impurities  in  the  commercial 
grade  driver  gas  in  the  region  where  the  driver  and  test 
gases  diffuse, together.  When  the  shock  luminosity  is  greater 
this  spike  vanishes  by  comparison. 

The  large  hump  (from  C  to  D)  is  due  to  scattered  light 
from  the  very  bright  plasma  behind  the  reflected  (n  =  2) 
shock.  This  light  grows  as  the  reflected  slug  lengthens  and 
approaches  window  #1,  then  decreases  again  as  the  interface 
causes  the  slug  to  decrease  in  size.  The  maximum  light 
here  Is  only  about  0.003  of  that  actually  seen  in  the  re¬ 
flected  shock  plasma  (second  picture).  As  shock  strength 
is  Increased  the  size  of  the  hump  (due  to  reflected  shock 
luminosity)  decreases  relative  to  the  maximum  primary  lumin¬ 
osity,  but  the  hump  is  always  present. 

For  stronger  shocks,  is  shorter,  but  the  larger 
shock  velocities  cause  the  hump  to  occur  much  earlier,  even 
before  the  interface  passes.  This  hump  on  top  of  the  flat 
primary  luminosity  curve  eventually  makes  it  impossible 
to  be  sure  that  relaxation  is  complete^  and  the  shock 
strength  range  to  be  studied  is  effectively  limited.  Treat¬ 
ing  the  test  section  walls  to  cut  the  reflections  would  help, 
but  treatment  with  molten  sodium  dichromate,  for  example, 
has  been  observed  to  cause  boundary  layer  problems  .  A 
better  solution  would  be  a  greater  separation  of  the  end- 
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wall  and  window  #1  so  that  the  scattered  light  from  the  re¬ 
flected  shock  will  be  reduced  considerably,  and  occur  much 
later. 

The  second  picture  shows  the  reflected  (n  =  2)  shock  at 
window  #2  (at  85  >asec.),  with  the  intensity  reduced  by 
neutral  density  filters.  The  filtering  has  cut  the  primary 
luminosity  and  tube  light  reflections  to  insignificant 
values  (X10  ).  Note  that  there  is  also  an  n  =  3  shock 

(at  135  nsec.),  apparently  reflected  from  the  interface 
just  ahead  of  window  #2,  since  the  interface  passes  at  150 
nsec.  Both  of  these  shocks  show  radiation  decay  (luminosity 
decrease  with  time). 

For  mercury,  the  luminosity  of  very  strong  shocks  be¬ 
came  so  large  that  observations  with  the  basic  system  could 
not  be  made.  The  luminosity  was  reduced  by  neutral  density 
filters  so  the  photomultiplier  was  not  overdriven  but  this 
reduced  the  schlieren  light  at  the  same  time.  When  the 
schlieren  light  is  reduced  too  far,  the  pulse  cannot  be  seen 

and  the  oscilloscopes  do  not  trigger  at  "zero"  time.  This 

* 

problem  was  overcome,  however,  by  using  a  laser  as  the 
schlieren  light  source  for  this  part  of  the  work.  The  laser 
is  extremely  bright,  which  allows  a  greater  luminosity 
range  to  be  studied  with  the  aid  of  neutral  density  filters. 
Color  filters  may  also  be  introduced  to  cut  the  mercury 
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luminosity  without  affecting  the  size  of  the  schlieren  pulse. 
The  schlieren  pulse  size  varies  with  p# ,  which  prohibits  the 
study  of  very  low  pe  shocks --unless  the  more  sensitive  laser 
schlieren  system  is  used. 


4.  Relaxation  Times 

The  laboratory  relaxation  time  (  )  is  not  the  relax¬ 

ation  time  (7")  which  the  gas  actually  experiences.  In  fact, 
the  time  scale  of  the  entire  luminosity  profile  is  distorted 
because  we  observe  gas  flowing  past  a  fixed  position,  in¬ 
stead  of  moving  with  the  gas  and  observing  the  relaxation. 
Zero  time  for  the  laboratory  time  scale  is  determined  by  the 
passage  of  the  shock.  At  laboratory  time  t* =  £/V,  (32) 

where  V  is  the  shock  velocity,  the  shock  will  be  Ic m  beyond 
the  window.  The  gas  which  is  then  at  the  window  was  actually 
shocked  t  seconds  earlier,  and  has  receded  from  the  shock  at 
a  relative  velocity  of  V-U,  where  U  is  the  gas  flow  velocity. 
Therefore,  the  particle  time  corresponding  to  tt  above  is 

t  =  i/(V-U).  (33) 

Eliminating  i  between  (32)  and  (33)  gives 

t  =  Vtk  /( V-U )  .  (3^) 


Using  Eq.  1  from  Ch.  II,  we  can  write  the  relation 
between  the  laboratory  and  particle  times  for  any  point  in 
the  relaxation  as 


t  = 


and, 


(  /po  )  1 1 

x,£  . 


* 


(35) 
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Note  that  this  derivation  assumes  that  U  (or  /),  )  remains 
constant  during  relaxation — this  is  only  approximately  true. 
The  laboratory  time  scale  is  not  only  incorrect  by  a  multi¬ 
plicative  factor,  but  is  actually  somewhat  distorted  due  to 
the  variations  of  at  (or  U). 


The  laboratory  relaxation  time  will  be  the  time  be¬ 
tween  the  schlieren  pulse  and  the  onset  of  constant  lumin¬ 
osity,  however,  determination  of  the  exact  time  at  which 
constant  luminosity  begins  is  usually  quite  difficult.  All 
oscilloscope  traces  show  some  electronic  noise  and  random 
fluctuations,  and  near  maximum  luminosity  the  relaxation 
curve  approaches  its  final  height  somewhat  asymptotically 
(see  M-10  Sec.  IV-2).  If  the  final  approach  is  truly  asymp¬ 
totic  then  theoretically  all  relaxation  times  are  infinite  — 
In  practice,  relaxation  is  complete  when  the  difference 
between  the  luminosity  and  the  final  asymptotic  height  is 
less  than  the  average  plasma  (and  noise)  fluctuation. 

In  the  present  research,  the  laboratory  relaxation  time 
is  taken  to  be  the  time  from  the  schlieren  pulse  until  the 
luminosity  Is  within  5$  of  its  final  asymptotic  value.  The 
final  height  can  be  determined  easily  and  quite  accurately 
for  all  pictures.  With  this  particular  choice  for  £  we 
are  not  usually  within  the  average  plasma  and  noise  fluct¬ 
uation  level  but  do  include  almost  all  of  the  relaxation 
curve  in  our  measurement.  Values  for  2)  >  which  are  ob- 


53 


talned  in  this  manner,  are  quite  consistent  when  a  picture 
is  measured  several  times  by  different  persons. 

5.  Errors  and  Corrections 

The  oscilloscope  sweeps  were  calibrated  using  a  1  MHz 
crystal  oscillator.  At  the  same  time  a  correction  was 
established  for  the  non-linearity  of  each  oscilloscope 
sweep.  All  sweep  corrections  were  quite  small  (<1%),  but 
were  applied  to  all  data  collected. 

The  use  of  three  windows  allows  shock  velocity  deter¬ 
mination  from  windows  #1  to  #2,  and  from  #2  to  #3,  to 
measure  the  shock  deceleration.  No  deceleration  was  noted 
for  any  argon  shock  but  all  mercury  shocks  showed  some, 
with  deceleration  becoming  greater  as  the  shock  strength 
increased.  The  effect  was  never  more  than  2$,  and  a  linear 
correction  was  applied  to  all  mercury  shocks  in  order  to 
find  the  shock  velocity  at  the  position  of  window  #1.  Cal¬ 
culations  allowed  for  the  thermal  expansion  of  the  tube. 

A  traveling  microscope  permitted  measurement  of  the 
pictures  to  better  than  -  0,01  mm.  With  such  measuring 
accuracy  the  shock  speeds  were  determined  to  within  i'  0.05$. 
The  tube  temperature  was  determined  to  better  than  t  0.1$. 
The  greatest  error  results  from  the  initial  pressure  values 
which  could  only  be  measured  slightly  better  than  t  1$ 

(for  p,  =  10  Torr--the  greatest  error). 

Errors  in  the  laboratory  variables  (U,  Pe  ,  T0 )  result 
in  errors  in  all  calculated  thermodynamic  variables  for  the 
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final  gas  conditions  (P;  ,  T, ,  etc).  Test  computations  show 
that  in  the  present  research  all  calculated  pressures,  mass 
densities,  and  particle  number  densities  will  be  within 
~  2$  of  the  correct  values  (error  due  mainly  to  4p#). 

Pinal  temperatures  are  in  doubt  by  t  3/4$  (due  to  errors  in 
T„  and  U).  The  values  for  electron  number  density  (Ne), 
electron  pressure  (Pe),  and  the  degree  of  ionization  (  ) 

are  the  most  sensitive  to  changes  in  the  initial  values,  but 
will  still  be  correct  to  within  t  4$. 

Laboratory  values  for  the  relaxation  times  can  always  be 
determined  within  *  3$j  with  most  within  t  2$--  the  accur¬ 
acy  depending  upon  the  smoothness  and  size  of  the  luminosity 
profile. 


CHAPTER  V 


EXPERIMENTAL  RESULTS 

A.  Mercury 
1.  Introduction 

The  data  from  the  present  study  have  been  reduced,  and 
the  results  are  displayed  in  a  series  of  graphs.  The  orig¬ 
inal  laboratory  data  (Pc,  T0  ,  V,  £)  are  tabulated  in  the 
Appendix.  The  shocks  observed  cover  a  range  of  Mach  numbers 
from  4.7  to  10,  and  Pfl  from  10.7  Torr  to  41.5  Torr.  Pinal 

plasma  temperatures  range  from  4035  °K  to  9626  °K,  with  , 

-5 

tie  final  degree  of  ionization  between  1.1  X  10  and  0.12. 
The  particle  relaxation  times  observed  are  from  10.8  jusec. 
to  548  psec. 

6,10 

It  is  common  practice  to  assume  that  the  relaxation 
process  will  obey  a  relationship  similar  to  many  chemical 
reactions,  that  is  ^  =  A  exp(E*/T),  (36) 

where  E*  is  the  "activation"  or  characteristic  energy  for 
the  process  (in  °K),  and  A  may  be  a  function  of  any  (or  all) 

of  the  plasma  thermodynamic  variables.  If  the  plot  of 
log  (G^)  against  T-1,  where  G  is  some  function  of  the 
thermodynamic  variables,  shows  less  point  scatter  than  the 
plot  of  log  T  against  T~\  then  T'  must  be  a  function  of  G. 

Therefore,  T"  =  B/G  exp(E/T),  (37) 

where  B  is  a  slower  varying  function  of  the  plasma  variables 
than  A.  If  G  includes  all  of  the  functional  relationships, 
then  B  will  be  a  constant,  and  the  graph  log(G^)  against 
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T  will  show  only  experimental  scatter.  The  slope  of  the 
graph  x 2. 303  will  be  E*  in  degrees  Kelvin. 

In  order  to  more  clearly  demonstrate  the  expected  de¬ 
pendence  of  T  on  P# ,  the  graphs  presented  distinguish 
shocks  on  the  basis  of  their  P«  values  (high,  medium,  low). 

The  curves  for  the  graphs  are  drawn  by  eye  to  represent  the 
best  fit  to  the  experimental  points.  Slope  variation  limits 
are  determined  by  the  most  extreme  curves  which  can  be  drawn 
consistent  with  the  data  points.  All  starred  variables 

■g- 

(eg.  T  )  refer  to  the  "ideal"  plasma  conditions,  that  is, 
the  conditions  Immediately  behind  the  shock  after  relaxation 
of  only  the  translational  states.  Unstarred  variables  (eg. 

T)  refer  to  the  final  (relaxed)  plasma  conditions.  All 
"Ideal"  and  final  variables  are  computed  as  indicated  in 
Chapter  II. 

2.  Pressure  Dependence  of  T 

Figures  10  and  11  correspond  to  the  graphs  of  Petschek 
and  Byron10  for  argon.  Figure  10  is  the  plot  of  log  ^  against 

l/T*.  iphg  plot  shows  a  tendency  for  points  to  form  a  pattern 

* 

sloping  up  to  the  right  (low  T  )  but  the  scatter  of  points 
is  very  extensive  (100-200$).  However,  from  the  identifica¬ 
tion  of  data  points  by  their  P0  value,  it  is  apparent  that 
high  P0  points  all  fall  low  in  the  plot  pattern,  while  low 
P0  points  are  all  high  In  the  pattern. 

# 

Figure  11  Is  the  plot  of  lot,  P,£  against  l/T  .  This 
plot  shows  considerably  less  scatter  (now  H0-70>)  than 
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10.  Plot  of  log?  against  1/T*  for  Mercury 
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Figure  11. 


Plot  of  log  ?t+  against  l/T*  for  Mercury 


Figure  10,  with  points  representing  different  P®  values 
fairly  well  intermingled  along  a  line  with  positive  slope. 
The  remaining  point  scatter  is  still  far  greater  than  that 
due  to  experimental  errors  (see  Sec.  1V-5)  hut  the  scatter 
has  been  reduced  significantly  by  the  introduction  of  P© . 

We  can  therefore  conclude  that  mercury  shows  the  same  de¬ 
pendence  on  initial  pressure  as  other  monatomic  gases — 
data  are  consistent  with  ^  =  BP^  exp(E*/T),  where  B  is 

Independent  of  P© ,  only  when  n  is  approximately  -1.  This 

6 

is  a  direct  contradiction  of  Kim's  result  . 

At  first  it  may  seem  surprising  that  T*  depends  upon  P«  , 
the  pressure  in  the  tube  before  the  shock  arrives.  For  the 

JL  # 

same  T0  and  T  ,  P©  is  proportional  to  N  ,  the  ideal  atom 

number  density  behind  the  shock.  Therefore,  Figure  11 

_  *  27 

actually  shows  that  C  depends  on  N  ,  and  some  researchers 

prefer  to  refer  to  this  as  a  density  dependence. 

3.  Temperature  Dependence  of  ^ 

5 

For  the  same  energy,  electrons  are  approximately  10 
times  as  effective  as  mercury  atoms  in  causing  ionization 
during  a  collision,  therefore  when  the  degree  of  ionization 
(©<)  becomes  greater  than  10~5  the  remaining  relaxation 
process  will  be  dominated  by  electron-atom  collisions. 
Petschek  and  Byron10  were  the  fJ-st  to  divide  the  total 
relaxation  time  into  stages  during  which  different  processes 
govern  the  relaxation.  They  were  also  the  first  to  point 
cut  that  the  electron  temperature  must  be  considerably  below 
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the  atom  temperature  during  relaxation. 

Different  temperatures  for  atoms  and  electrons  in  the 

same  gas  are  possible  because  the  two  components  are  well 

insulated  from  one  another  energetically.  The  fraction  of 

a  heavy  atom's  energy  which  is  transferred  to  an  electron 

in  a  collision  is  roughly  equal  to  the  mass  ratio  of  the 

particles  For  the  electron -mercury  atom  collisions  this 

-6 

is  appx-oximately  10”  .  Since  the  average  atom  thermal 

energy  is eV,  or  about  one-tenth  of  the  ionization 
_7 

potential,  10  of  the  energy  required  to  create  a  new 

electron  is  transferred  to  the  electrons  in  each  collision. 

7 

Therefore,  10  collisions  are  necessary  for  each  new 
electron  produced,  but  only  one  or  two  collisions  between 
electrons  will  cause  a  Maxwellian  distribution.  The  elect¬ 
rons  will  therefore  be  in  Maxwellian  distribution  at  all 
times,  with  the  high  energy  tail  of  the  distribution  con¬ 
tinually  losing  energy  as  ions  are  created. 

Harwell  and  Jahn2®  later  showed  experimentally  that  for 
argon,  xenon,  or  krypton,  the  early  stage  of  ionizational  re¬ 
laxation  was  due  mainly  to  atom-atom  collisions.  Atom- 
atom  collisions  first  create  exc  ted  atoms,  then  further 
collisions  of  the  excited  atoms  result  in  ionization.  Ion¬ 
ization  is  controlled  by  the  production  rate  of  the  excited 
atoms.  Many  people  have  attempted,  with  limited  success, 
to  describe  the  relaxation  process  mainly  in  terms  of  these 
two  types  of  collision  processes.  The  most  complete 
mathematical  treatment  Is  that  of  Wong  and  Bershader  ,  who 
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also  calculate  the  electron  temperature  (Te),  relative  to 
the  atom  temperature  (T* )  at  all  times  during  relaxation. 

30 

Prom  their  work,  and  the  experiments  of  Morgan  and  Morrison  , 
we  know  that  the  electron  temperature  is  approximately 
equal  (but  slightly  lower)  to  the  final  plasma  temperature 
(T)  throughout  relaxation. 

For  relaxations  where  the  final  degree  of  ionization  is 

high  (as  in  the  present  research),  the  electron-atom  collision 

process  will  be  dominant  throughout  much  of  the  relaxation 

* 

time.  Therefore  Te  (=■  T),  and  not  T  should  be  the  relevant 

* 

temperature  for  the  relaxation.  The  use  of  T  would  only 
be  correct  for  studies  of  the  initial  ionization  rates2^’^0 
or  where  only  the  atom-atom  collision  process  occurs. 

Figure  12  is  the  plot  of  log  T  vs  T  1  which  shows  essent¬ 
ially  the  same  features  as  Figure  10.  Comparing  Figures  10 

* 

and  12  it  can  be  seen  that  the  use  of  T  rather  than  T 
causes  the  pressure  dependence  to  appear  less  significant 
at  the  .lower  (high  T)  end  of  the  graph --even  though  the 
range  of  P  is  greatest  here.  Figure  13  is  the  plot  of 
log  PoT  vs  T_1  which  shows  a  fairly  good  line  with  positive 
slope  over  most  of  the  range  but  with  more  scatter  and  in¬ 
dication  of  a  bend  (towards  the  horizontal)  at  the  low  T  end 
of  the  plot.  Figure  13  shows  significantly  less  scatter 
than  the  corresponding  Figure  11  especially  at  the  high 
temperature  end  of  the  graph. 

From  the  comparison  of  Figures  10  and  11  with  Figures 
12  and  13  we  can  conclude  that  the  final  plasma  temperature 
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(or  T c)  is  a  more  significant  temperature  than  T*  for  the 
total  relaxation  time  in  a  highly  ionized  plasma. 

4.  Radiation  Decay 

Both  Figures  11  and  13  seem  to  show  an  over-correction 
for  pressure  in  the  lowest  part  of  the  graph  (low  P0  values 
are  below  the  medium  %  values).  This  is  probably  due 
to  the  strong  radiation  decay  which  is  present  behind  the 
high  temperature  shocks.  Shocks  typical  of  this  high  temp¬ 
erature  region  are  shown  in  Figure  14.  The  attainment  of 
equilibrium  is  taken  to  be  the  highest  point  of  the  trace, 
although  the  form  of  the  curve  suggests  that  if  radiation 
decay  were  not  present,  the  luminosity  might  have  risen 
higher  than  this  value.  Therefore,  the  measured  values 
for  Tl  at  high  temperatures  may  be  too  short — the  error 
increasing  as  the  effect  of  the  radiation  decay  increases 
with  temperature.  Pictures  from  shocks  M54  and  M55  show 
the  increase  of  radiation  decay  with  final  temperature. 

Shock  m66,  at  a  higher  temperature,  actually  shows  a  decay 
between  the  otner  two  shocks.  It  would  appear  that,  for  the 
same  final  temperature,  radiation  decay  is  more  important 
for  lower  pressures.  This  is  probably  due  to  the  increased 
optical  thickness  (especially  to  resonance  radiation)  of 
the  higher  pressure  gas.  We  might  thus  expect  that  the 
measured  for  low  P0  (high  T)  points  in  Figures  11  and 
13  are  decreased  more  than  the  higher  P„  (high  T)  points. 


This  radiation  effect  would  also  tend  to  account  for  the 
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apparent  lack  of  pressure  dependence  which  Kim  observed. 

All  of  his  shocks  fall  in  the  very  high  temperature  and  low 
P o  region  where  the  effect  is  greatest  (8000-10,000  °K 
and  lo-20  Torr). 

5.  Dependence  of  Z"  on  Electron  Number  Density 

The  comparison  of  the  relaxation  process  with  chemical 
reactions  which  resulted  in  Eq.  37  was  quite  superficial 
and  should  be  re-examined.  For  the  type  of  chemical  re¬ 
action  which  most  nearly  resembles  the  collision-ionization 

31 

process  Fowler  and  Guggenheim  show  that  the  reaction  rate 
constant  (number  of  reactions  per  cc  per  sec.)  is, 

k  =  C'1  exp(-E„/T)  ,  (38) 

where  E„  is  the  activation  energy,  and  C  may  be  a  function 
of  the  chemical  concentrations  and/or  the  temperature  (the 
plasma  variables  in  our  case). 

In  the  present  reaction  k  =  dNe/dt  ,  where  Ne  is 
the  electron  number  density  (or  atomic  ion  density)  at  any 
Instant.  We  may  therefore  re-write  Eq.  38  and  integrate 
over  the  relaxation  process  to  give, 

A  A 

J dt  =  /  C  exp{E/i/T)dNe  ,  (39) 

<5  Q 

where  Ne  is  the  final  electron  number  density  in  the  plasma. 
Eq.  39  can  be  integrated  to  give, 

Z*  =  CNe  exp(E„/T),  (**0) 

if  T  and  c  remain  constant  throughout  the  relaxation--which 
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is  approximately  true293}  (from  before  Te  =  T  final). 

Prom  Eq.  40  we  would  expect  to  be  a  function  of  Ne. 

Figure  15  is  the  plot  of  log(f/Ne)  against  T-1.  In 
Figure  15  the  pressure  dependence  is  quite  apparent.  The 
plot  separates  distinctly  into  points  grouped  on  the  basis 
of  their  P0  values,  with  high  P0  shocks  always  lowest.  The 
point  scatter  may  have  been  reduced  by  Ne  (re  Fig.  12),  but 
this  is  not  immediately  apparent  because  of  the  different 
graph  scale  used. 

Figure  16  is  the  plot  of  log(N*?-/Nt)  against  T"1.  This 
graph  shows  an  excellent  grouping  of  points,  independent  of 
their  P0  values,  along  a  curve  with  positive  slope.  The 
curve  appears  to  be  a  straight  line  which  bends  toward  the 
horizontal  at  lower  temperatures.  A  careful  comparison  of 
Figures  13  and  16  shows  that  the  introduction  of  Ne  has  re¬ 
duced  the  point  scatter  veAy  slightly  (by~10$)  especially 
at  the  low  T  end  of  the  graphs.  ;re  may  thus  conclude  that 
does  depend  upon  the  electron  number  density.  Since 
electron-atom  collisions  play  a  major  role  In  the  present 
relaxations  it  is  not  surprising  to  find  that  't'  depends  on 
both  N*  and  N  . 

Figure  17,  log(N*7'/Ne)  against  (T*),  is  plotted  for 
comparison  with  Figure  16,  to  demonstrate  that  with  the  Ne 
dependence  of  T  included,  the  use  of  T  rather  than  T*  still 
results  in  less  scatter. 


10 


Figure  15.  Plot  of  log(?/Ne)  against  1/T  for  Mercury 


Figure  16.  Plot  of  log(N*?VNe)  against  1/T  for  Mercury. 
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6.  Activation  Energy 

From  Figure  16,  we  may  now  say  that 

T  = C(Nc/N*)  exp(Ew/T)  , 

where  C  is  a  constant,  and  E*  apparently  varies  in  the 
temperature  range  covered  (note  change  in  slope).  The 
activation  energy  should  be  the  energy  of  che  rate -controlling 
step  of  the  relaxation  reaction,  but  two  collision  reactions 
occur  during  relaxation  which  may  have  different  values  of 
E/,  .  For  the  high  temperature  end  of  Figure  16,  where  re¬ 
laxation  depends  mainly  on  electron-atom  collisions,  the 
slope  of  the  graph  gives  E*  =  (11.5  t  1.0)  x  102*  °K, 

which  corresponds  quite  well  with  the  ionization  potential 
of  the  mercury  atom  (^  =  12.1  x  10^  °K  or  10.43  eV). 

This  result  indicates  that  the  electron-atom-dominated 
relaxation  proceeds  in  a  single  step  rather  than  by  inter- 
mediate  stages  as  Kim  concluded  . 

The  variation  in  E*  (slope)  doe3  not  necessarily  mean 
that  E*  is  temperature  dependent.  The  value  obtained  will 
be  the  average,  over  the  entire  *£■  ,  of  the  E*  effective  at 
each  instant  during  relaxation.  For  the  atom-atom  collision 
phase  of  lonizational  relaxation,  E*  has  been  found  to 
correspond  to  .tie  first  excited  state  of  the  atom  for  argon, 
krypton,  xenon 2®  and  cesium2^.  Thus,  E*  from  Figure  16 
should  be  slightly  low,  because  the  first  part  of  the  re¬ 
laxation  must  proceed  by  atom-atom  collisions,  which  depend 
upon  the  first  excited  state  of  the  atom  (5.42  x  10^  °K  or 
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4.76  eV).  For  lower  temperature  shocks  (upper  end  of  Fig.  16) 
the  atom-atom  process  will  occupy  an  increasing  fraction  of 
the  relaxation,  so  we  expect,  and  find,  even  lower  values 
of  l^(*=7.3  x  102*  °K  at  the  top). 

If  the  Ne  dependence  is  not  accepted,  then  from  Figure  13 
using  only  the  low  ?0  points,  we  may  obtain  a  value  for 
E*(  =  5.2  x  104  °K)  approximately  equal  to  the  energy  of  the 
first  excited  state  of  the  atom,  ill  other  values  of  E* 
determined  from  Figures  11  or  13  are  too  low--E*  must  be 
equal  to,  or  greater  than,  the  energy  of  the  first  excited 
state  of  the  mercury  atom.  Values  of  E„  less  than  the  energy 
of  the  the  first  excited  state  cannot  be  explained  by  im¬ 
purity  effects,  since  the  excited  states  of  most  impurities 
will  be  above  that  of  the  mercur’y,  and  the  amount  of  impurity 
is  far  too  low,. 

B.  Argon 

Figures  13  and  19  show  the  Importer''  plots  for  argon. 

again  depends  on  P0  (N*)  and  N» ,  and  the  scatter  <~f  ex¬ 
perimental  points  was  observed  to  be  less  for  T  than  for 
T  ,  as  with  mercury.  From  Figure  19,  E*  =  (11.4  _  O.o) 
x  10^  °K  whicn  is  less  tnan  the  energy  of  the  first  excited 
state  of  the  argon  atom  (13.4  x  lC2i),  but  approximately 
equal  to  the  Ionization  potential  of  mercury.  There  may 
have  been  apDreciable  mercury  vapor  present  as  an  impurity 
in  all  argon  shocks.  It  was  not  realized  at  the  tJme,  but 


Figure  18.  Plot  of  3og  P/f  against  l/T  for  Argon 


Figure  19.  Plot  of  log(N*?/Ne)  against  1/T  for  Argon 


if  the  water  cooled  trap  (between  the  shock  tube  and  the 
diffusion  pump)  was  not  efficient  enough,  the  mercury  dif¬ 
fusion  pump  might  actually  maintain  several  microns  of 
mercury  vapor  in  the  heated  tube  at  all  times.  Therefore, 
as  pointed  out  by  Jones  and  McChesney  ,  the  observed  E* 
should  be  that  for  an  easily  excited  impurity,  even  for 
quite  low  impurity  concentrations.  The  fairly  high  mercury 
contamination  will  act  as  the  source  of  electrons  for  the 
argon  reaction  and  the  mercury  ionization  energy  will  be 
seen  as  the  effective  E„ . 

Because  of  the  suspected  mercury  contamination  of  all 
argon  shocks,  conclusions  based  solely  on  the  present  argon 
data  are  not  considered  reliable. 


CHAPTER  VI 


SPECIAL  RESULTS  RELATED  TO  THE  LUMINOSITY  PROFILE 

1.  Radiation  Decay 

Pictures  of  luminosity  traces  for  three  strong  mercury 
shocks  are  shown  in  Figure  20.  The  shape  of  these  curves, 

7 

agrees  extremely  well  with  those  of  Tumakaev  and  Lazovskaya  . 
Their  curves  are  actually  plots  which  represent  the  measured 
populations  of  individual  mercury  atom  excited  states  at 
various  distances  behind  the  shock  The  observed  total 
luminosity  profile  apparently  agrees  with  that  due  to  indi¬ 
vidual  spectral  lines — at  least  for  strong  shocks. 

Tumakaev  and  Lazovskaya  also  claim  that  a  considerable 
number  of  excited  atoms,  due  to  photoexcitation  resulting 
from  diffusion  of  resonance  radiation,  are  present  in  the 
mercury  vapor  ahead  of  the  shock  wave.  In  the  present  study, 
very  strong  luminosity  was  detected  before  the  arrival  of 
the  shock,  but  whether  this  was  radiation  from  excited  atoms, 
or  light  scattered  from  the  walls  of  the  tube,  is  not  known. 
This  pre-luminosity  was  so  great  for  strong  shocks  that  the 
original  schlieren  light  source  was  too  faint  to  be  detected 
and  the  brighter  laser  schlieren  system  was  needed.  Weymann 

and  Holmes  have  studied  similar  pre -luminosity  due  to  elect¬ 
ro  3^ 

ron  diffusion-5-3,  and  due  to  photoionization  ,  in  argon 
where  the  effect  is  considerably  smaller  than  In  mercury. 

In  Figure  20  the  upper  beam  trace  of  M5^  clearly  shows  that 
the  luminosity  at  window  #2  is  not  zero  when  the  shock 
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Figure  20.  Radiation  Decay  Curves  for  Mercury. 
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arrives. 

6,25 

Kim  states  that  boundary  layers  are  very  thick  for 
all  strong  mercury  shocks.  This  would  suggest  that  a  large 
part  of  the  apparent  radiation  decay  may  actually  be  the 
result  of  observing  the  plasma  through  growing  boundary  lay¬ 
ers.  Laminar  boundary  flows  are  difficult  to  detect.  With 
the  sensitive  laser  schlieren  system,  however,  some  traces 
show  extra  spikes  (  eg.  M65)  which  can  only  be  from  the 
schlieren  system.  This  indicates  that  turbulence  is  present, 
either  in  the  plasma  or  in  a  boundary  layer. 

During  relaxation,  radiation  decay  may  be  relatively 
unimportant  until  the  very  end  (and  after),  but  boundary 
layers,  which  absorb  material  from  the  plasma,  could  have 
a  very  strong  influence  on  the  observed  relaxation  time. 

The  question  of  whether  the  obseived  luminosity  decay  is  due 
to  boundary  layers  or  radiation  decay,  could  be  answered  by 
observing  the  same  strong  shock  simultaneously  through  diff¬ 
erent  plasma  thicknesses — this  is  possible  with  a  rectangular 
tube  cross-section.  The  radiation  decay  rate  would  be  un¬ 
affected,  while  the  ratio  of  plasma  to  boundary  layer  thick¬ 
ness  would  definitely  change. 

2.  Unusual  Features  of  the  Luminosity  Profile 

The  traces  of  the  three  different  shocks  given  in  Figure 
21  demonstrate  a  rather  unexpected  luminosity  characteristic 
which  is  present  for  both  argon  and  mercury.  Each  profile 
contains  two  distinct  points  at  which  the  slope  of  the  curve 
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changes  abruptly.  These  bends  a  n  be  seen  in  alL  shock 
traces  (except  very  strong  mercury  shocks)  but  are  often  not 
as  pronounced  as  in  the  sample  tr- ces.  No  reason  could  be 
found  for  why  the  bends  are  more  -ronoui.ced  in  some  shocks 
than  in  others  and  the  reason  for  this  unusual  luminosity 
feature  is  not  known.  A  better  understanding  of  this  feature 
and  why  it  occurs  may  aid  in  understanding  the  entire  relax¬ 
ation  process. 

This  feature  of  relaxation  has  not  been  noted  previously, 
possibly  because  no  previous  research  used  this  particular 
method  of  observing  the  luminosity.  The  total  luminosity 
is  usually  recorded  directly  on  film,  with  a  different 
spectral  range  being  recorded  (usually  much  further  into  the 
U.V. ),  and  much  poorer  time  (space)  resolution  than  the 
present  system.  The  bends  have  apparently  not  been  observed 
for  the  curves  of  growth  of  the  electron  density  during  re¬ 
laxation,  using  microwave  reflection  or  electron  probe  tech¬ 
niques.  Neither  of  these  techniques,  however,  can  observe 
the  entire  relaxation  curve  efficiently,  and  again  the  time 
resolution  is  not  generally  as  good  as  the  present  system. 

Thus  we  may  either  assume  that  this  feature  is  actual¬ 
ly  present  but  was  obscured  by  poor  resolution  in  other  ex¬ 
periments,  or,  it  is  the  result  of  observing  a  limited  and 
integrated  spectral  range  with  the  present  technique. 

The  second  alternative  does  not  seem  likely  since  the 
total  radiation  should  increase  uniformly  during  relaxation 
and  there  is  no  apparent  reason  why  one  spectral  range  (or 


Figure  22.  Plots  of  Normalized  Luminosity  Profiles 
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strong  line?)  should  vary  differently  fr^m  the  total.  This 
question  could  be  decided  by  an  experiment  which  would  ob¬ 
serve  the  total  luminosity  and  a  single  spectral  line,  say 
the  resonance  line  2537  A°,  simultaneously  for  comparison. 
Soma  possible  experimental  causes  for  the  unusual  curve 

QC. 

shape  can  be  eliminated.  Pox,  McLaren  and  HobsonJV  have 
given  the  transformation  from  laboratory  to  particle  time, 
and  their  work  eliminates  the  laboratory  time  scale  dis¬ 
tortions  (much  too  small)  as  the  source  of  the  unusual 
shape.  Several  shock  traces,  obtained  without  the  presence 
of  any  schlieren  system  still  show  bends  (see  Pig.  22), 
thus  eliminating  the  schlieren  system  as  the  cause  of  this 
phenomenon . 

There  are  some  indications  In  other  research  results 

that  this  feature  is  actually  present  in  all  shocks. 

Kfcj-ly’s^  Figure  5  does  seem  to  show  the  same  effect  for 

microwave  measurements  of  the  electron  density  in  an  argon 

29 

shock.  The  experimental  points  of  Wong  and  Bershader 
vary  from  their  theoretical  curves  in  a  manner  also  slightly 

suggestive  of  this  same  feature.  Petschek  and  Byron10 
briefly  observed  the  argon  continuum  using  a  phototube,  as 

verification  for  their  electron  probe  technique,  and  although 
their  curves  show  a  great  deal  of  noise  fluctuation  there 
is  a  suggestion  of  the  same  bends. 

Figure  22  shows  the  plots  for  several  traces  which  have 
been  measured  and  normalized  (  7" =  1  and  maximum  height  =  1) 
for  comparison.  From  these  graphs,  and  consideration  of  the 
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other  traces,  it  appears  that  the  second  ^end  always  occurs 
between  48$  and  55$  of  the  total  relaxation  time  (7l). 
Therefore,  the  time  to  the  second  bend  (to)  satisfies 
(though  not  as  well)  ail  the  same  relationships  as  T  (Pigs. 

13  and  1 6).  No  reasonable  relationship  could  be  found  be¬ 
tween  ti,  the  time  until  the  f,  ’St  bend  occurs,  and  the 
plasma  variables — perhaps  because  the  measurements  of  t-^ 
are  not  very  accurate  (±  30$). 

The  time  from  to  to  ^  may  correspond  to  the  electron- 
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atom  phase  of  relaxation--Haught  indicates  that  this  phase 
occupies  approximately  IT-*  Perhaps  ti  is  the  time  for  the 
atom-atom  phase,  vith  t-^  to  t2  representing  tne  region 
where  both  collision  processes  are  competing,  or  an  unknown 
third  process  exists. 

3.  Impurity  Effects 

Figure  23  shows  the  effect  of  small  amounts  of  air  im¬ 
purity  (deliberately  Introduced  via  tube  leak)  on  the  mer¬ 
cury  luminosity  profile.  For  extremely  low  impurity  levels 
the  characteristic  "overshoot"  form  appears.  The  luminosity 
overshoots  the  final  equilibrium  value,  and  then  returns  to 
equilibrium  shortly  later,  thus  producing  a  hump  at  the  be¬ 
ginning  of  the  constant  luminosity  region.  In  the  present 
research,  tests  prove  tnat  radiation  overshoot  is  definitely 
a  low  impurity  level  effect.  For  impurity  levels  below  that 
of  M24  (0.005$)  no  profile  deformation  occurs.  The  same  form 
appears  in  argon  but  at  an  impurity  level  about  four  times 
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Figure  23.  Radiation  Overshoot 
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higher  than  for  mercury.  Haught2?  claimed  very  high  gas 
purity  for  his  study  of  cesium  (-<1  part  per  million),  but 
he  observed  radiation  overshoot  for  all  of  his  shocks.  It 
would  appear  that  impurities  are  more  important  when  the 
first  excited  state  of  the  atom  is  lower. 

Morgan  and  Morrison  have  criticized  the  use  of  the 
luminosity  delay  time  as  the  relaxation  time,  mainly  on  the 
basis  of  Haught’s  luminosity  curves.  However,  the  problem 
is  only  to  decide  at  which  point  on  the  trace  the  relaxation 
is  actually  complete.  When  the  length  of  the  overshoot  hump 
is  appreciable  relative  to  Ti  ,  the  best  point  to  use  is 
easily  determined.  From  the  traces  for  several  shocks  it 
has  been  determined  that  if  TL  is  taken  as  the  beginning  of 
the  flat  region  Immediately  behind  the  hump  then  the  shock 
fits  extremely  well  on  Figures  13  and  16.  This  is  the  same 
point  which  Haught  arbitrarily  chose  for  his  work.  T  is 
still  deflined  as  the  time  until  the  luminosity  is  within 
5$  of  its  final  asymptotic  value. 

Figure  24  shows  three  shock  traces  from  a  series  taken 
to  demonstrate  the  effect  of  impurities  on  argon--all  shocks 
were  approximately  the  same,  except  for  the  amount  of  im¬ 
purity.  Actually,  two  series  were  run,  one  with  hydrogen 
as  the  impurity  and  the  other  with  air,  but  the  two  series 
are  indistinguishable.  AI-2  has  the  lowest  impurity  level 
w.ilch  shows  an  effect--the  first  hump  is  overshoot,  the 
second  is  due  to  scattered  light  from  the  reflected  shock. 
For  higher  impurity  levels  the  overshoot  hump  becomes  less 
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apparent  as  the  entire  profile  becomes  greatly  distorted. 
Impurities  affect  the  relaxation  zone  only,  although  the 
entire  trace  seems  to  become  more  erratic  (noisy)  when  im¬ 
purities  are  present. 

Figure  25  shows  some  types  of  contamination  which  could 
accidently  occur.  MI-3  demonstrates  dibutyl  pthalate 
(D.B.P. )  vapor  contamination.  D.B.P.  is  used  in  the  manomet¬ 
er  for  measuring  gas  pressures.  In  the  present  research, 
the  D.B.P.  vapor  is  isolated  from  the  shock  tube  by  the 
bellows  gauge,  but  if  gas  pressures  are  measured  directly 
using  an  oil  manometer,  contamination  will  result.  This  is 
extremely  important  because  many  researchers  use  oil  manomet¬ 
ers  for  pressure  measurements  and  then  claim  unusually  high 
gas  purity,  whereas  oil  contamination  may  actually  be  in¬ 
fluencing  all  of  their  results. 

Small  amounts  of  impurities  are  often  deliberately  added 
to  the  test  gas  in  order  to  identify  the  shock  position  by 
a  large  luminous  (molecular)  spike.  MI-3  shows  that  the  oil 
molecular  spike  is  actually  about  4  ^seconds  behind  the 
shock  (for  a  strong  shock).  An  error  of  this  size  may  only  be 
ignored  in  experiments  where  T  is  very  large.  The  actual 
time  error  for  a  given  impurity  can  only  be  determined  by 
experiment,  and  will  obviously  depend  upon  the  type  of  im¬ 
purity  used. 

AI-10  shows  Duo-seal  pump  oil  contamination.  At  the 
slower  sweep  speed  no  delay  is  noted  for  the  molecular  spike. 


Figure  25.  Oil  Contamination  Profiles 


but  the  spike  at  zero  time  is  much  too  broad  to  te  the 
schlieren  system  alone.  This  particular  oil  corf aninat ion 
was  introduced  accidently.  When  the  oil  in  th  *ao-seal 
pump  for  the  driver  section  became  too  low,  the-  '"ii  overheat¬ 
ed  and  oil  vapor  was  introduced  into  the  driver  section  be¬ 
tween  shocks.  The  oil  contamination  was  then  carried  into 
the  expansion  chamber  when  a  shock  was  fired,  and  the  next 
shock  showed  oil  contamination. 

Noting  the  contamination  associated  wit.i  two  types  of 
oil,  an  attempt  was  made  to  introduce  contamination  due  to 
a  common  diffusion  pump  oil  (Convoil-20  Vacuum  Pump  Fluid). 
Many  shock  tubes  use  oil  diffusion  pumps  instead  of  mercury. 
Under  no  circumstances  could  any  significant  affect  on  the 
luminosity  profile  be  produced. 

The  mercury  contamination  of  argon  chocks  did  noticeably 
affect  the  luminosity  profile  but  did  change  the  total  re¬ 
laxation  time  significantly.  The  same  result  could  also  be 
true  of  any  type  of  impurity  contamination.  Mercury  impurity 
in  argon  did  not  cause  overshoot--perhaps  because  the  ion¬ 
ization  potential  of  mercury  is  slightly  lower  than  the  first 
excitation  level  of  argon,  rather  than  being  higher  as  in 
the  case  of  nitrogen  or  hydrogen  impurity. 


CHAPTER  VII 


SUMMARY 


Comparison  of  the  point  scatter,  in  the  plots  of  log  ^ 
against  l/T*  and  log  against  1/T*  for  mercury,  estab¬ 

lishes  that  ionizational  relaxation  time  depends  approx¬ 
imately  on  Pc  1.  The  fact  that  log  (P,T  )  against  1/T 
results  in  even  less  scatter  demonstrates  that  for  highly 

ionized  relaxations,  the  final  plasma  temperature  T  is  more 

* 

relevant  than  the  commonly  used  ideal  gas  temperature  T  . 
Point  scatter  is  further  reduced  very  slightly  in  the  plot 
of  log( P„  ^  /Ne )  against  1/T  indicating  that  ionizational 
relaxation  time  is  approximately  proportional  to  Ne,  the 
final  electron  density  of  the  plasma.  The  Indus  Jon  of  Ne 
also  results  In  a  more  acceptable  value  for  the  activation 
energy. 

The  value  obtained  for  the  activation  energy  indicates 
that  Ionization  resulting  from  electron-atom  collisions  pro¬ 
ceeds  as  a  single  step  rather  than  by  cumulative  steps  as  in 
the  case  of  atom-atom  collisions. 

The  above  conclusions  may  also  be  drawn  from  the  argon 
data  but  the  value  obtained  for  the  activation  energy  in¬ 
dicates  that  all  argon  shocks  may  contain  appreciable  mer¬ 
cury  as  an  impurity.  For  this  reason  conclusions  are  not 
based  solely  on  the  argon  data. 

It  is  not  likely  that  further  research,  which  compares 
the  total  ionizational  relaxation  times  of  different  shocks. 


can  contribute  significantly  to  a  better  understanding  of 
the  relaxation  process.  Future  studies  of  ionizatlonal 
relaxation  should  attempt  to  observe  individual  shocks  in 
greater  detail  rather  than  attempting  to  relate  many  diff¬ 
erent  shocks.  The  ionizational  relaxation  phenomenon 
appears  to  be  too  complex  to  be  analysed  by  any  one  obser¬ 
vation  technique  alone.  Each  individual  shock  should  be  ex¬ 
amined  simultaneously  using  several  different  techniques  and 
the  results  related  if  the  relaxation  processes  are  ever  to 
be  completely  understood.  There  is  more  information  contain¬ 
ed  in  each  shock  than  can  be  obtained  from  a  single  number 
representing  the  total  relaxation  time. 

The  observed  luminosity  profiles  of  the  present  research 
show  two  bends.  This  is  the  first  time  this  feature  has  been 
observed,  probably  because  of  the  unusually  high  spatial 
resolution  of  the  present  observation  technique..  (0.006  inch). 
The  very  poor  spatial  resolution  of  drum  camera  techniques 
(0.150  inch)  makes  them  virtually  useless  for  studying  the 
relaxation  process.  The  reason  for  the  bends  in  the  lumin¬ 
osity  profile  is  not  known.  The  shape  of  the  profile  sug¬ 
gests  that  the  total  luminosity  consists  of  at  least  two 
components  which  vary  differently  in  form.  It  is  suggested 
that  a  future  study  should  compare  individual  spectral  lines 
with  one  another,  with  the  continuum  and  with  limited 
spectral  regions,  using  the  present  observation  techniques. 
Special  attention  should  be  given  to  the  resonance  line 
2537  A0.  It  would  be  helpful  to  the  understanding  of  re- 
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laxation  processes  to  know  how  different  spectral  regions 
and  lines  compare  during  the  relaxation. 

For  strong  shocks  in  mercury,  radiation  decay  is  appar¬ 
ently  very  pronounced;  however,  there  are  indications  that 
part  of  the  observed  effect  is  actually  due  to  boundary 
layers.  The  proportion  of  the  effect  which  is  due  to  bound¬ 
ary  layers  can  be  determined  by  comparing  the  same  shock 
simultanecusly  through  different  plasma  thicknesses  in  a 
rectangular  cross-section  tube.  If  boundary  layers  are  very 
pronounced  then  the  results  of  studies  in  the  strong  mer¬ 
cury  shock  region  will  always  be  questionable. 

Radiation  overshoot  is  a  low  impurity  level  effect  and 
the  relaxation  time,  for  shocks  showing  overshoot,  appears 
to  correspond  to  the  end  of  the  overshoot  hump.  Air  and 
hydrogen  as  impurities  appear  to  greatly  affect  the  shape 
of  the  relaxation  luminosity  profile. 

Manometer  oil,  or  pump  oil  contamination  can  produce 
very  pronounced  changes  in  the  relaxation  processes.  This 
is  extremely  important  because  it  shows  the  need  to  guard 
against  these  contaminations  in  all  future  experiments  and 
also  casts  doubt  on  the  results  of  some  previous  experiments. 


APPENDIX 


The  measured  laboratory  variables  for  the  relaxation 
studies  of  this  research  are  given  for  future  reference. 

The  mercury  data  are  tabulated  in  Table  I,  and  although  the 
argon  shocks  are  believed  to  be  contaminated  by  mercury  vapor 
the  argon  data  are  tabulated  in  Table  II. 


The  following  laboratory  data  values  are  given: 

V  The  shock  velocity  at  window  #1,  in  cm  per  psecond. 

Corrections  have  been  made  for  the  thermal  expansion  of 
the  tube,  and  for  deceleration  in  the  mercury  shocks. 

?0  The  initial  pressure  in  she  tube  ahead  of  the  shock, 
in  Torr,  was  measures  in  mm  of  D.B. P.  but  values  have 
been  converted  to  Torr  using  the  measured  density  of  D.B. P. 

T0  The  temperature  of  the  tube  ahead  of  che  shock.  In  °C, 

'k  The  laboratory  relaxation  time-  .n  pseeonds,  Is  taken  as 

Lu  the  time  from  shock  passage  until  the  luminosity  is 
within  5$  of  its  final  asymptotic  value. 

t:  The  laboratory  time  from  shock  passage  until  the  first 
bend  in  the  luminosity  profile,  in  p&econds. 

to  The  laboratory  time  from  shock  passage  until  the  second 
bend  in  the  luminosity  profile.  In  pseconds. 


Corrections  have  been  made  to  V,  £  >  ancJ  ^2  -0I> 

oscilloscope  sweep  rate  differences  and  non-linearity. 
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Table  1.  Mercury  Laboratory  Dcta 


SbocK 

V 

( cm/sec ) 

Po  % 
(Torr) 

To 

(°c) 

(nsec) 

ti 

(}isec) 

t2  x 
(nsec) 

Ml 

M2 

M3 

M4 

111,090 

103,420 

115,680 

111,750 

16.57 

18.19 

l6.ll 

26.05 

242 

242 

242 

241 

36.8 

86.7 

30.7 
20.0 

13  0 

3.8 

44.3 
15.7 

11.3 

m6 

M7 

m8 

M9 

116,600 

108,700 

111,090 

100,900 

10.71 

19.34 

16.03 

16.11 

238 

237 

2^2.5 

240 

38.9 

46.8 

40.6 

94.2 

4.9 

9.5 

18.2 

27.7 

23.5 

Mio 

Mil 

Ml  2 

Ml  3 

100,820 

37,940 

104,420 

92,920 

16,11 

14.18 

12.49 

20.96 

240 

241 

241 

241 

95.9 

135.1 

148.8 

110.4 

10.7 

26.4 

19.9 

46.3 
67.1 

69.3 
59.0 

Ml  4 

M15 

Ml6 

M17 

92,470 

107,290 

99,120 

106,140 

18.65 

13.87 

16.49 

15.72 

241 

240.5 

240 

240 

114.5 

81.5 

104.1 

79.2 

10.8 

6.4 

10.7 

58.8 

38.7 

50.6 

40.1 

mi8 

M19 

M20 

M21 

101,210 

109,760 

99,850 

100,440 

18.03 

15.49 

14.72 

14.87 

244 

244 

245 
244.5 

97.2 
54.6 

86.2 
103.8 

7.4 
11.7 

7.6 

7.5 

45.4 

30.2 

38.7 

55.9 

M22 

M2  3 
M24 
M25 

101.980 

110,160 

112,930 

109,070 

14.57 

14.26 

15.67 

15.34 

245.5 
2^5 

245 

244.5 

98.6 

67.0 

34.6 

48.4 

12.7 

6.4 

8.5 

53.7 

35.2 

16.0 

24.9 

M26 

M28 

M29 

M30 

98,660 

89,050 

91,220 

90,500 

14.95 

30.83 

30.83 

30.60 

245 

244 

244 

244 

126.9 

90.4 

61.8 

60.7 

10.  ( 

6~4 

60,  7 
47.3 

32.1 

M31 
M3  3 
M35 
M36 

95,070 

104,940 

300,130 

105,830 

33.60 

33.29 

34.83 

12.41 

244 

252 

250 

246 

53.1 

34.2 

44.3 
75.6 

7.5 

5.3 

10.7 

5.8 

28.1 

15.7 

38.5 

M38 

M39 

M4l 

M42 

101,050 

97,130 

111,600 

103,020 

32.37 

34.14 

29.0.5 

31.06 

2M3 

242.5 

241 

242 

46.5 

43.4 

13.05 

49.2 

7.5 

22.9 

6.8 

24.1 
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Table  1.  Mercury  Laboratory  Data  (con't). 


Shock 

V 

Po 

T0 

7i 

ti 

( cm/sec ) 

(Torr) 

(°c) 

(Msec) 

(Msec) 

(usee) 

M43 

105,630 

29.29 

242 

28.6 

m45 

94,850 

15.88 

2l  2 

137.1 

_ 

59.9 

M46 

100,880 

17.26 

242 

93.7 

60.1 

M47 

101,970 

15.41 

242 

85.2 

6.4 

40.8 

M48 

107,000 

15.26 

240 

83.3 

4.4 

43.2 

M49 

101,910 

16.65 

240 

71.2 

6.4 

38.3 

M50 

107,150 

16.18 

240 

60.8 

M51 

107,480 

15.57 

240 

64.3 

11.0 

35.8 

M52 

108,740 

16.03 

242 

53.7 

6.6 

28.1 

M53 

126,290 

15.41 

242 

20.5 

12.0 

M54 

138,370 

12.33 

242 

10.57 

4.8 

M55 

145,540 

11.71 

242 

8.49 

- 

4.3 

M56 

160 , 270 

12.87 

242 

4.16 

2.0 

M57 

173,390 

11.87 

242 

2.85 

M58 

173,730 

12.72 

242 

2.75 

M59 

129,770 

30.98 

243 

6.43 

- 

- 

m6o 

131,090 

32.98 

243 

5.33 

m6i 

120,680 

31.14 

243 

9.28 

_ 

m62 

120,690 

32.52 

243 

9.87 

_ 

M63 

140,350 

30.52 

243 

3.76 

- 

- 

m64 

141,530 

41.00 

242 

2.23 

M65 

140,410 

37.15 

242 

2.4s 

m66 

144,630 

41.46 

240 

2.12 

_ 

__ 

M67 

161,660 

22.20 

240 

2.75 

- 

- 

m68 

172,890 

18.34 

238 

2.65 

Mo9 

187,880 

17.57 

240 

1.925 

— 

— 

M70 

156,130 

15.57 

240 

3.54 

- 

— 

M7l 

144,490 

16.03 

240 

5.78 

- 

- 

U)  ro 


96 


Table  2.  Argon  Laboratory  Data 


Shock 

V 

Po 

To 

ti 

(cm/sec) 

(Torr) 

(°c) 

(usee) 

(nsec) 

(nsec) 

Al 

296,580 

19.22 

198 

44.8 

6.3 

23.0 

A2 

340,470 

9.43 

198 

39.9 

9.0 

21.0 

A3 

293,680 

13.65 

198 

49.1 

8.9 

24.2 

A4 

320,550 

21.42 

198 

44.2 

9.1 

18.9 

A5 

315,520 

21.12 

195 

42.0 

9.9 

23.2 

a6 

305,440 

21.42 

195 

38.5 

5.6 

21.5 

A7 

282,400 

21.50 

200 

54.8 

12.7 

23.4 

a8 

236,700 

23.35 

197 

52.7 

- 

A9 

299,820 

19.53 

197 

47.4 

7.5 

26.6 

A10 

309,040 

19.73 

197 

47.0 

5.6 

25.0 

All 

313,060 

21.12 

197 

46.4 

8.4 

24.1 

A12 

318,850 

19.88 

197 

45.9 

8.8 

21.6 

A13 

327,520 

19.34 

198 

39.0 

11.4 

19.3 

A14 

332,030 

18.80 

198 

35.5 

11.4 

18.7 

A15 

336,570 

18.22 

163 

43.4 

21.0 

A16 

337,510 

18.80 

165 

45.0 

11.2 

20.6 

A17 

304,270 

27.74 

164 

41.0 

11.9 

22.2 

ai8 

319,610 

26.20 

165 

46.4 

10.9 

21.8 

A20 

338,460 

14.80 

167 

'35.3 

11.0 

20.7 

A21 

3s 4, 260 

13.26 

167 

34.1 

- 

A22 

313,780 

27.51 

166 

42.7 

11.3 

21.8 

A24 

305,820 

26.51 

162 

28.1 

5.8 

16.0 

A25 

322,690 

19.27 

163 

29.2 

2.3 

13.2 

A26 

297,660 

16.65 

163 

46.2 

6.0 

21.8 

A27 

276,990 

19.88 

163 

51.7 

— 

28.1 

A29 

237,940 

17.73 

163 

65.5 

- 

32.4 

A  30 

240,030 

19.19 

162.5 

51.1 

- 

27.3 

A31 

259,710 

17.65 

162 

62.0 

- 

A32 

251,590 

18.80 

160 

142.8 

45.4 

81.5 

A33 

261,350 

19.96 

160 

71.3 

- 

- 

A34 

258,840 

21.19 

162 

57.3 

- 

25.5 

A  37 

270,340 

20.11 

161 

249.0 

- 

113.5 

A  39 

286,210 

29.36 

162 

50.3 

9.9 

24.8 

A40 

292,460 

29.13 

161 

40.0 

11.2 

19.0 

A4l 

279,560 

32.29 

160 

46.7 

6.7 

21.3 

A42 

275,100 

39.30 

160.5 

50.0 

6.3 

28.0 

§ 

t 


97 


Table  2.  Argon  Laboratory  Data  (con't). 


Shock 

V 

Po 

To 

Tl 

ti 

tg 

( cm/sec) 

(Torr) 

(°c) 

(jasec) 

(»sec) 

(usee) 

A44 

285,260 

21.58 

162 

45.0 

7.0 

24.5 

A46 

322,170 

21.89 

162.5 

35.5 

8.6 

A47 

298,990 

32.29 

161.5 

42.7 

9.2 

23.8 

A54 

301,490 

40,69 

254 

45.6 

10.5 

23.1 

A56 

331.310 

26.67 

254 

46.3 

25.1 

A57 

315,700 

29.05 

255 

41.3 

6.2 

21.2 

A58 

299,390 

46.09 

254 

40.2 

5.5 

21.0 

A59 

249,680 

38.15 

252.5 

54.1 

3.9 

23.8 
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